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n-Particle Irreducible Effective Actions (nPIEA) are a powerful tool for extracting non- 
perturbative and non-equilibrium physics from quantum field theories. Unfortunately, practical 
truncations of nPIEA can unphysically violate symmetries. Pilaftsis and Teresi (PT) addressed this 
by introducing a “symmetry improvement” scheme in the context of the 2PIEA for an O (2) scalar 
theory, ensuring that the Goldstone boson is massless in the broken symmetry phase [A. Pilaftsis 
and D. Teresi, Nuclear Physics B 874, 2 (2013), pp. 594-619.]. We extend this idea by introducing 
a symmetry improved 3PIEA for O (N) theories, for which the basic variables are the one-, two- 
and three-point correlation functions. This requires the imposition of a Ward identity involving the 
three-point function. We find that the method leads to an infinity of physically distinct schemes, 
though a field theoretic analogue of d’Alembert’s principle is used to single out a unique scheme. 
The standard equivalence hierarchy of nPIEA no longer holds with symmetry improvement and we 
investigate the difference between the symmetry improved 3PIEA and 2PIEA. We present renor¬ 
malized equations of motion and counter-terms for two and three loop truncations of the effective 
action, though we leave their numerical solution to future work. We solve the Hartree-Fock approx¬ 
imation and find that our method achieves a middle ground between the unimproved 2PIEA and 
PT methods. The phase transition predicted by our method is weakly first order and the Goldstone 
theorem is satisfied, while the PT method correctly predicts a second order phase transition. In 
contrast, the unimproved 2PIEA predicts a strong first order transition with large violations of the 
Goldstone theorem. We also show that, in contrast to PT, the two loop truncation of the symmetry 
improved 3PIEA does not predict the correct Higgs decay rate although the three loop truncation 
does, at least to leading order. These results suggest that symmetry improvement should not be 
applied to nPIEA truncated to < n loops. We also show that symmetry improvement schemes are 
compatible with the Coleman-Mermin-Wagner theorem, giving a check on the consistency of the 
formalism. 
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I. INTRODUCTION 

The recent demands of non-equilibrium field theory ap¬ 
plications in particle physics, cosmology and condensed 
matter have led to a renaissance in the development of 
novel field theory methods. The ^-matrix school, re¬ 
booted in the guise of spinor-helicity methods, has led to 
a dramatic speedup in the computation of gauge the¬ 
ory scattering amplitudes in vacuum [I|. On the fi¬ 
nite temperature and density fronts, efficient functional 
methods in the form of n-particle irreducible effective 
actions (nPIEA) have proven useful to understand col¬ 
lective behaviour and phase transitions 0. They are 
similar in spirit to methods based on Schwinger-Dyson 
equations in field theory or BBGKY (Bogoliubov-Born- 
Green-Kirkwood-Yvon) equations in kinetic theory how¬ 
ever, unlike the Schwinger-Dyson or BBGKY equations, 
nPIEA naturally form closed systems of equations of mo¬ 
tion without requiring any closure ansatz [H3. nPIEA 
methods can be understood as a hybrid of variational 
and perturbative methods: nPIEA consist of a series of 
Feynman diagrams, however the propagators and ver¬ 
tices of these diagrams are the exact 1- through n-point 
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proper connected correlation functions which are deter¬ 
mined self-consistently using variational equations of mo¬ 
tion. 

This self-consistency effectively resums certain classes 
of perturbative Feynman diagrams to infinite order. For 
example, the one loop 2PIEA diagram corresponding to 
the Hartree-Fock self-energy in theory actually sums 
all of the so-called daisy and super-daisy graphs of or¬ 
dinary perturbation theory (Figure ll.l|) . This particular 
resummation is often done in the literature without the 
use of nPIEA, but such ad hoc resummation schemes run 
the risk of summing an asymptotic series: a mathemat¬ 
ically dangerous operation (recent progress on summa- 
bility has been made in resurgence theory [^, which is 
beyond the scope of this work). nPIEA sidestep this issue 
because they are defined by the rigorous Legendre trans¬ 
form procedure, guaranteeing equivalence with the orig¬ 
inal theory. Unlike ad hoc resummations, nPIEA based 
approximation schemes are placed on a firm theoretical 
footing and can be systematically improved. 

However, loop-wise truncations of nPIEA, n > 1, 
have difficulties in the treatment of theories with sponta¬ 
neously broken continuous symmetries. The root cause 
of these difficulties is the fact that nPIEA obey different 
Ward identities than the IPIEA. When the effective ac¬ 
tion is truncated to a finite order the equivalence between 
the Ward identities is lost. This can also be understood 
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FIG. I.l. From left to right: the Hartree-Fock self-energy di¬ 
agram, an example daisy or ring diagram, an example super¬ 
daisy graph. The whole class of super-daisy diagrams is ob¬ 
tained from iterating insertions of Hartree-Fock graphs in all 
possible ways. 


in terms of the resummation of perturbative Feynman di¬ 
agrams: when an nPIEA is truncated some subset of per¬ 
turbative diagrams are summed to infinite order, but the 
complementary subset is left out entirely. The pattern 
of resummations does not guarantee that the cancella¬ 
tions between perturbative diagrams needed to maintain 
the symmetry are kept. In the case of scalar field theo¬ 
ries with O {N) —7> O {N — 1) breaking, the result is that 
the final O (iV — 1) symmetry is maintained, but, at the 
Hartree-Fock level of approximation, the nonlinearly re¬ 
alised O (N) /O {N — 1) is lost, the Goldstone theorem is 
violated (the — 1 Goldstone bosons are massive), and 
the symmetry restoration phase transition is first order in 
contradiction with the second order transition expected 
on the basis of universality arguments. A similar prob¬ 
lem arises in gauge theories, where the violation of gauge 
invariance in the Gloop truncation is due to the missing 
{I + l)-loop diagrams (see, e.g. 0, 0 for a discussion of 
the gauge fixing problem). 

Several studies have attempted to find a remedy for 
this problem. These are discussed in @ and references 
therein. Here we restrict attention to the technique most 
frequently advocated in the literature [a. This tech¬ 
nique constructs the so called external propagator as the 
second functional derivative of a resummed effective ac¬ 
tion which depends only on the mean field, obtained by 
eliminating the 2- through n-point correlation functions 
of the nPIEA by their equations of motion. The result¬ 
ing effective action does obey a IPI type Ward identity 
and the external propagator yields massless Goldstone 
bosons. However, the external propagator is not the 
propagator used in loop graphs, so the loop corrections 
still contain massive Goldstone bosons leading to incor¬ 
rect thresholds, decay rates and violations of unitarity. In 
order to avoid these problems a manifestly self-consistent 
scheme must be used. 

Pilaftsis and Teresi recently developed a method which 
circumvents these difficulties @ for the widely used 
2PIEA (also known as the CJT effective action after 
Gornwall, Jackiw and Tomboulis [HI, the Luttinger- 
Ward functional or <I>-derivable approximation depending 
on the context). The idea is incredibly simple: impose 
the desired Ward identities directly on the free correla¬ 
tion functions. This is consistently implemented by us¬ 
ing Lagrange multipliers. The remarkable point is that 
the resulting equations of motion can be put into a form 
that completely eliminates the Lagrange multiplier field. 


They achieve this by taking a limit in which the Lagrange 
multiplier vanishes from all but one of the equations of 
motion, and this remaining equation of motion is replaced 
with the constraint to obtain a closed system. We show 
that this nontrivial aspect of the procedure generalizes 
to the 3PIEA. We find that the generalization requires a 
careful consideration of the variational procedure, how¬ 
ever, and an infinity of schemes are possible. A new 
principle is required to choose between the schemes and 
we propose what we call the d’Alembert formalism as 
the appropriate principle by analogy to the constrained 
variational problem in mechanics. 

We extend the work of Pilaftsis and Teresi to the 
3PIEA for three reasons. First, the 3PIEA is known to 
be the required starting point to obtain a self-consistent 
non-equilibrium kinetic theory of gauge theories. The ac¬ 
curate calculation of transport coefficients and thermal- 
ization times in gauge theories requires the use of nPIEA 
with n > 3 (see, e.g. d, [13, US] and references therein 
for discussion). The fundamental reason for this is that 
the 3PIEA includes medium induced effects on the three- 
point vertex at leading order. The 2PIEA in gauge theory 
contains a dressed propagator but not a dressed vertex, 
leading not only to an inconsistency of the resulting ki¬ 
netic equation but also to a spurious gauge dependence of 
the kind discussed previously. We consider this work to 
be a stepping stone towards a fully self-consistent, non- 
perturbative and manifestly gauge invariant treatment of 
out of equilibrium gauge theories. 

Second, nPIEA allows one to accurately describe the 
initial value problem with 1- to n-point connected cor¬ 
relation functions in the initial state. For example, the 
widely used 2PIEA allows one to solve the initial value 
problem for initial states with a Gaussian density matrix. 
However, the physical applications one has in mind typ¬ 
ically start from a near thermal equilibrium state which 
is not well approximated by a Gaussian density matrix. 
This leads to problems with renormalization, unphysi¬ 
cal transient resp onses and thermalization to the wrong 
temperature [13]. This is addressed in [T3| by the addi¬ 
tion of an infinite set of nonlocal vertices which only have 
support at the initial time. Going to n > 2 allows one 
to better describe the initial state, thereby reducing the 
need for additional nonlocal vertices. 

Lastly, the infinite hierarchy of nPIEA is the natural 
home for the widely used 2PIEA (in all its guises) and 
provides the clearest route for systematic improvements 
over existing treatments. Thus investigating symmetry 
improvement of 3PIEA is a well motived next step in the 
development of non-perturbative QFT. 

After this introductory section we review nPIEA in 
Section m focusing on the 3PIEA for a model O (N) 
scalar field theory with symmetry breaking as a specific 
example. Then in Section m we review and extend the 
symmetry improvement program of Pilaftsis and Teresi. 
This includes a derivation of the required Ward identi¬ 
ties, their implementation as constraints using Lagrange 
multipliers and the limiting procedure required to ob- 
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tain sensible equations of motion for the system. We 
will see that this procedure rests on a certain techni¬ 
cal assumption which we will justify in Appendix lAl and 
make a connection to the d’Alembert principle using a 
mechanical analogy. Then in Section IIVI we investigate 
the renormalization of the theory, first with the two loop 
truncation and then three loops. The three loop trun¬ 
cation is analytically intractable in 1 -|- 3 dimensions so, 
after discussing the renormalization procedure in arbi¬ 
trary dimension, we present results for 1-1-2 dimensions. 
The result of this section is a set of finite equations of 
motion which must be solved numerically. In Section |V] 
we solve the theory at the Hartree-Fock level and discuss 
the phase transition thermodynamics. Section I VII is a 
verification that the Coleman-Mermin-Wagner theorem 
holds in the symmetry improvement formalism despite 
the imposition of Ward identities, a check on the con¬ 
sistency of the formalism. In Section IVIII we discuss the 
effects of symmetry improvement on the absorptive parts 
of propagators and make some comments involving the 
Higgs decay rate and dispersion relations. Finally in Sec¬ 
tion IVIIII we discuss the main themes of the paper and 
point out directions for future work. 


II. REVIEW OF nPI EFFECTIVE ACTIONS 

For the sake of having an explicit example we consider 
the O (N) linear cr-model given by the action 

S [(!>]= J - ^rri^4ia4>'' - ^ , (n.l) 

where a = 1, • • • , N is the flavor index. In the symme¬ 
try breaking regime rri^ < 0 and a vacuum expectation 
value develops, which by symmetry can be taken in the 
last component {(j)) = (0,...,0,n) where = —Qm^lX 
at tree level. The massive mode, which we loosely call 
“the Higgs” (reflecting our ultimate interest in the Stan¬ 
dard Model, despite the absence of gauge interactions in 
(111.11) 1. gets a tree level mass = An^/3 = — 2 to^. 

The nPI effective actions form a systematic hierarchy 
of functionals F^") [</?, A, V, • • • , where (/?,••• > 
are the proper 1- through n-point correlation functions 
and we have suppressed spacetime arguments and flavour 
indices. In more detail, 

= i^a) , (H.2) 

Aab =ih{{T [(t)a(t>b]) - {(pa) {(pb)) , (H.3) 

h?IXad^be^cfVdef = (T [(pa(pb(pc]) “ (T [(pa(pb]) {(pc) 

- (T [(pc(pa]) {(pb) - (T [(pb(pc\) {(pa) 

+ 2 {bPa) {(Pb) {(Pc) (H.4) 


On notation: we work mostly in 1-1-3 dimensions with 
Vd-c ~ diag (1, — 1, — 1, — 1), although the generalization 
to other dimensions is simple. We take h = c = i 
as far as units are concerned, though we keep loop 
counting factors of h explicit. Repeated indices are 
summed. Often, field indices accompany spacetime ar¬ 
guments. Repeated indices in this case imply an inte¬ 
gration over the corresponding spacetime argument as 
well (“DeWitt notation”). Where explicitly indicated, 
spacetime and momentum integrals are written in com¬ 
pressed notation with f^ = f d^x, fp = f d'^p/ (27r)^ and 

/p = /p d^p/ (27r)^ etc. (T [• • • ]) represents the time or¬ 
dered product of the factors in [•••]. Through most of 
this article the meaning of time ordering is left implicit. 
The formalism can be readily applied to vacuum field 
theory {t G (—oo,-|-oo) with the natural ordering), fi¬ 
nite temperature field theory in the imaginary time or 
Matsubara formalism {t —>■ —ir, with periodic boundary 


conditions on r G 


0, P = and the natural ordering 

on r) 1^, and general non-equilibrium field theory on 
the two time Schwinger-Keldysh contour {t runs from 0 
to -|-oo then from -|-oo — ie back down to 0 — ie with time 
ordering in the sense of position along the contour rather 
than the magnitude |t|) [l^. In Section |TV] we develop 
the renormalization theory for the vacuum case and in 
Section s we solve the Hartree-Fock approximation at 
finite temperature in the Matsubara formalism. 


In general 1/^”) is the sum of connected one particle irre¬ 
ducible Feynman diagrams contributing to (0") with all 
external legs (including leg corrections) removed. 

In the absence of external source terms the correlation 
functions obey equations of motion of the form 




= 0 , 




= 0 , 




= 0 . 


(H.5) 


6(p ’ (5A ’ ’ i5I/l" 

In the exact theory F*)”) obey equivalence relationships 

r(i) = r(2) A] = r(3) [p, a, v] = • • •, (ii.e) 


where extra arguments are eliminated by their equations 
of motion when comparisons are made. These relation¬ 
ships only hold approximately when approximations are 
made to the theory. A stronger equivalence hierarchy 
that relates loop-wise truncations of the F^”) will be dis¬ 
cussed below. 

For later convenience we introduce the tree level vertex 
functions 


Voabc {x,y,z) = 
Wabcd {x,y,z,w) = 


6^S[(P] 


5(pa {x) S(pb (y) S(pc{z) 

S^S[^] 


(H.7) 


p—if 


Ha (x) 5(pb {y) He {z) Hd {w) 


(H.8) 

















4 


For the O (N) model these are 

Voabc {x, y,z) = - ^ [Sabfc {x) + Sca^b (x) + Sbc^^a (a;)] 

X (a; — y) {x — z), (II-9) 

hFabcd ix^ y^ Z^ ra) — ~ {^ab^cd 4“ ^ac^bd 4“ ^ad^bc] 

X (x — y) {x — z) {x — w). 

( 11 . 10 ) 

I 


The nPIEA is defined in the functional integral for¬ 
malism by the Legendre transform of the connected gen¬ 
erating function 




n) 


( 11 . 11 ) 

for a held theory in the presence of source terms dehned 
by the generating functional 


J,iC(2),.. . ,X(") 


V [(j)] exp - ( 5 [0] -1- 


(j)y 


+ 

i I 

r?. 


( 11 . 12 ) 


Then r(’") is the n-fold Legendre transform 


r(") ip,A,v,-- -,- j- 


SJ 


-K^n) 


<5PF(”) 

Jk^ 

(5A:(") ' 


(11.13) 


where the source terms are solved for 

in terms of the (/?, A,-- - ,14(’^). Spacetime integrations 
and O (N) index contractions have been suppressed for 
brevity. For bosonic helds the A, • • • , are totally 
symmetric under permutations of their arguments. The 
generalisation to fermions requires sign changes for odd 
permutations of arguments corresponding to fermionic 
helds, but is otherwise straightforward. (Note that the 
nPIEA is dehned by this Legendre transform, not by any 
irreducibility property of the Feynman graphs, though for 
low enough loop orders the graphs are irreducible as the 
name implies. At high enough loop order for n > 2 the 
name becomes misleading. For example, the hve loop 
5PIEA contains graphs that are not hve-particle irre¬ 
ducible 0!) 

r(i) is the familiar IPI effective action introduced 
by Goldstone, Salam and Weinberg and independently 
by Jona-Lasinio [l^. F^^^ \if\ can be written 


F^i) M = 5 [(/s] -f yTrln {Aq ^45]}+ M , (11.14) 


where Ag ^ [tp] = 5'^S[(j) + ip] /is the inverse 

propagator and F^^^ is the sum of all connected vacuum 
graphs with > 2 loops where the propagators Aq [p] and 
vertices are obtained from the shifted action S [(j> + p] 
with the additional prescription that all 1-particle re¬ 
ducible graphs are dropped [T^ • Note that (jIL4ll is equiv¬ 
alent to 



53r(i) 

5pd (u) 6pe {v)Spf (w)’ 


(11.15) 


where the superscript “(1)” indicates the vertex derived 
from the IPIEA. 

The 2PIEA r(2) [(^, A] was introduced in the context 
of non-relativistic statistical mechanics, apparently in¬ 
dependently, by Lee and Yang, Luttinger and Ward, 
and others, but was brought to the functional formal¬ 
ism and relativistic field theory by Cornwall, Jackiw and 
Tomboulis [HI- F^^^ [p, A] is most easily computed by 
noting that the Legendre transform can be performed 
in stages. First perform the Legendre transform with 
respect to J, using result pi.141) with the replacement 
S'[(^] —>■ S'[(()] 4- \ct)xK^J(l)y, then do the transform with 
respect to . This procedure leads to (up to an irrel¬ 
evant constant) 

r(2) A] = S [p] + yTrln (A’l) + ^Tr (Ag-^A) 

+ f(")[(p,A]. ( 11 . 16 ) 

The equation of motion for A is Dyson’s equation. 


A-i = Ag-i - S, 

where 

^ 2i 5T^^\p,A] 

h 5A 


(11.17) 


(11.18) 


is identified as the IPI self-energy. Since E consists of 
( 2 ) 

IPI two-point graphs, F 2 must consist of 2PI vacuum 

(2i) 

graphs. That is, F 2 ' is the sum of all vacuum diagrams 
which do not fall apart when any two lines are cut. This 
results in a drastic reduction in the number of graphs at a 
given loop order. Further, the propagators in a 2PI graph 
are the full propagators A, with all self-energy insertions 
resummed to infinite order. 

The 3PIEA F^^^ [p^ A, V] can be computed following 

the same method: replace S' [i^] —>■ S' [4>] 4- ^K^yz4>x4>y4‘z 
in the previous result and perform the Legendre trans¬ 
form with respect to 
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FIG. II.1. Two and three loop diagrams contributing to 
Fg®^ We label these 3>i through $5 from left to 

right, respectively, and their explicit forms are given in (III.211) 
through (III.251) . Solid circles represent the resummed ver¬ 
tices V and the crossed circles represent the bare vertices 
Vo and W. The dashed lines represent the resummed prop¬ 
agators A. Note that these diagrams are called “eight,” 
“egg,” “mergedes,” “hair,” and “bball” respectively in 
the nomenclature of Bi- 


The shift by the source term K^^'> results in the 
introduction of an effective three point vertex V = 
Vq + K^^'> appearing in The difficult step of 

the Legendre transform is relating V to V. This can 
be done by comparing [j, with 


(5r(2) 
r(3) is 


ip,A;V 


/SK^^'> (see iEl). The final result for 


r(3) =S[^] + |Trln(A-i) + ^Tr (A^-iA) + T^^^, 

(11.19) 

where to three loop order the diagram piece is 

j-2 

rf^ = $1 -h ^VoAAAV - $2 

+ $3 + «>4 + $5 + o {n ‘^), ( 11 . 20 ) 


where - ,$5 are given by the Feynman diagrams 

shown in Figure UlTI Explicitly, 


permuted as well). The graphical interpretation of this 
equation is shown in Figure |lL2] The permutations lead 
to the usual s, t, and u channel contributions with the ex¬ 
pected symmetry factors. This equation is best thought 
of as a self-consistent integral equation in the same spirit 
as a Schwinger-Dyson equation, and can be solved it¬ 
eratively. By iterating (111.261) one sees that it sums a 
sequence of vertex correction diagrams to infinite order. 

If all higher order terms are kept in Fg the resulting V 
is the same as the of (III.151) . however truncated ac¬ 
tions give solutions V ^ Similar remarks apply for 
the propagators. These self-consistent solutions do not 
in general obey the desirable field theoretic properties of 
the full solution, such as Ward identities. The symmetry 
improvement strategy is to impose IPI Ward identities 
which as constraints on the self-consistent solutions A 
and V. This is discussed further in Section imi 

Note that pi. 261) can be derived by removing a re- 

/ oV 

summed vertex from each graph in Fg (because S/dV 
acts by removing a single V factor from graphs in all 
possible ways), which has the graphical effect of opening 
two loops. This means that the one loop correction to 
V comes from three loop graphs in F^^). Thus a loop- 
wise truncation of nPIEA for n > 3 does not lead to a 
loop-wise truncation of the corresponding equations of 
motion. We will discuss the further implications of this 
in Section I VIII 

Another important implication of this result is that 
F( 2 ) and F(3) are equivalent to two loop order (after one 
substitutes V = Vq + 0{h) in However, F^^^ and 

r(3) differ at three loop order because F^^) contains re¬ 
summed vertex corrections that F^^^ does not. This is 
an example of an equivalence hierarchy of nPI effective 
actions that has the general form 


‘hi = --^WabcdAabAcd, (11-21) 

O 

‘^2 = -r^VabcVdef AadAbeAcf, ( 11 . 22 ) 


‘hs — ^abc^def^ghi^jklAadAbgAcjAeliAfkAii^ 

(11.23) 

ih? 

‘1*4 — ^Vabci^def^^ghijAadAbgAQbAeiAfj^ (11.24) 
‘hs = -J^WabcdWefghAaeAbfAcgAdh- (11.25) 


The 3PI equation of motion for V is 0 = IfY"’ 
reads in full 

Pabc — boafec “f i^A^adei^hfg^chiAdfA^hAgi 

~ ^ 14-(o)delhV(b)7r(c)/g^d/^eg + O {k ) , 

(11.26) 


where is a sum over the 3! permutations mapping 
{a,b,c) —>■ (tt (a), TT (b), TT (c)) (spacetime arguments are 


r?ioop)['^] = rS\oop)[‘^>A] = ..., (11.27) 

rg\oop) [‘^] ^ rS 2 \oop) A] = [<(>, a, v] = - - -, 

(11.28) 


la^oop) [^] ^ A] ^ [^, A, V] 


k 2 ) 


(3) 


= r 


(4) 

(3 loop) 




(11.29) 


where the subscripts represent the order of the loop-wise 
truncation and the “extra” correlation functions are to 
be evaluated at the solutions of their respective equa¬ 
tions of motion before making the comparison (and also 
allowance is made for shifts by irrelevant constants). This 
equivalence hierarchy has been explicitly checked up to 
five loop 5PI order in scalar field theories Q. 

The existence of the equivalence hierarchy implies that 
in the standard formalism one gains nothing by going to 
higher nPI effective actions unless one also includes di¬ 
agrams with at least n loops, since for m > n one can 
always reduce to r(”\oop). However, we shall see 
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FIG. II.2. Equation of motion for the 3PI vertex function V 
up to one loop order (III.261) . Note that the bubble graph (last 
term) is implicitly symmetrized over external momenta and 
O (N) indices. 


that symmetry improvement breaks this equivalence hier¬ 
archy. In particular, we find that the symmetry improve¬ 
ment of the 3PI effective action modifies the A equation 
of motion in a way that remains non-trivial even if 
is then truncated at two loops and V is replaced by its 
tree level value Vq ■ In general we find that the symmetry 
improvement procedure introduces Ward identities that 
relate fc-point functions to (fc -I- l)-point functions and 
these constraints spoil the equivalence hierarchy, i.e. the 
“operations” of symmetry improvement and reduction in 


the hierarchy do not commute. The consequences of this 
for the phase diagram of the scalar O {N) theory in the 
various possible schemes are investigated in Section ©■ 

III. SYMMETRY IMPROVEMENT 

Symmetry improvement begins with the consideration 
of the Ward identities in the nPI formalism. Following 
Q we derive the Ward identities from the condition that 
the effective action is invariant under a symmetry trans¬ 
formation. The theory in (III. II) has the O (iV) symmetry 
transform 

0a </^’a + (IH-l) 

where are the generators of the group in the funda¬ 
mental representation (yl = 1, • • • ,N {N — 1) /2) and ca 
are infinitesimal transformation parameters. Note that 
our implicit integration convention can be maintained if 
we consider that contains a spacetime delta function 
oc (xq — Xb). Also = —Tf^. Under this trans¬ 
formation the effective actions change by 


= ^-^^eAT^b‘fb, 

^r(2) ^ ^ {T^Acb + niAac) , 

0(fa O/Xab 

ieAT^tpb + iCA {r^^Acb + T^Aac) + J^J—ieA iT^^Vdbc + T^Vadc + T^Vabd) , 


S(Pa 


ab 


6Vabc 


(111.2) 

(111.3) 

(111.4) 


according to the tensorial structure of the arguments. 
The next steps to derive the Ward identities are to set 
ijr(") = 0, take functional derivatives of the resulting 
equations with respect to ip and finally apply the equa¬ 
tions of motion. We also extract the overall factors of ieA- 
We call the identity derived from the m-th derivative of 
the (m + l)-point nPI Ward identity, denoted by 

yVa}—am = 0 where oi, • • • , Om are O (fV)/spacetime in¬ 
dices. We note first of all that = 0 identically by 
the equations of motion. We also find that 


y^A{i) 


jpd) 

SpcSpa 

sr(^ 


y'ab'-Pb, 


SpdSpcS'Pa 


y'ab^b + 




(111.5) 

(111. 6 ) 


Specialising now to the broken symmetry vacuum pb = 
v5bN we obtain the following identities by substituting 


different generators in turn: 


0=/ ^ca {Xc,Xa)v, N, 

0 = / VNab (x, y,z)v + A~d (x, y) 

J Z 

f^afeAjYjY (x, y), a,b ^ N 

0 = ^ca T bl C, 

0 = / I4ca (x, y, z) V, d,c,a^ N, 

J Z 

0 = / VNNa (x, y, z)v, N. 

J Z 


(III.7) 


(111. 8 ) 

(111.9) 

(111.10) 

(111.11) 


Note that we explicitly write spacetime arguments, 
O (N) indices and integrations in the above. This is be¬ 
cause DeWitt notation would lead to ambiguities here. 
Below we introduce an ansatz adapted to the situation 
which again allows for notational simplifications. 
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The essence of symmetry improvement is to impose 
these Ward identities, derived for the IPI correlation 
functions, on the nPI correlation functions. Effectively, 
we change / (Aipi, Eipi) / (A 3 P 1 , E 3 P 1 ), where / is the 
Ward identity and we change the arguments but not the 
functional form. We have already made this substitution 
in dnm - dnLTT]) . 

The first two identities will prove important in the fol¬ 
lowing, however, the identities (IIII.9I) - (IIII.11I) are trivial 
in the sense that they can be satisfied simply by postu¬ 
lating an ansatz for A and V which is tensorial under the 
unbroken O (JV — 1) symmetry. For later convenience we 
adopt this spontaneous symmetry breaking (SSB) ansatz 
now by introducing the notation 


r Ag (a:, 2/), a = b^ N, 

Aab{x,y) = IAh {x,y), a = h = N, (III.12) 
I 0 , otherwise, 

for the Goldstone (Ag) and Higgs {Ah) propagators re¬ 
spectively, and we also introduce the vertex functions V 


and Vn where 


Vabc {x, y, z) 


'V {x,y,z)6aNSbc 
-bcyclic permutations 

Vn (x, y, z) 

0 


exactly one of 
a,b,c = N, 

a = b = c = N, 
otherwise. 

(III.13) 


Note that Vn is not constrained by any of the 2 - or 
3-point Ward identities. Spacetime arguments are per¬ 
muted along with the O (N) indices, so that the first 
spacetime argument of V is always the one referring to 
the Higgs. The other two arguments refer to the Gold- 
stone bosons and V is symmetric under their interchange. 
Vn is totally symmetric in its arguments. For reference 
note that at the 2 -loop truncation, V = Vq and we obtain 
F = {-Xv/3) X (x - y) 5 V) {x - z) and Vn = 3H. Af¬ 
ter substituting the ansatz the diagrams $ 1 , • • • , $5 can 
be put into the form 


$1 = - 1 ) AgAg -b {N - 1) AgAh + —^AhAh, 

^2 = — {N — 1) VVAhAgAg + —VnVn {Ah)^ , 

$3 = (AT - 1) (H)^ (Ah)^ (Ag)^ + ^ {Vn)'^ {Auf + {N - 1) — (H)^ Ah Ah {Ag)^ , 


$4 


$5 


= 2{N — 1) VVn (Ah)^ Ag Ag -I- (A^ — l) VV Ah {Ag)‘^ + 3VnVn (Ah)^ 

+2" {N-1)VV{Ag)''AhAh] , 

= ^ {[(A - 1) AgAg + AnAnf + 2{N-1) (Ag)" + 2 (A^^)"} . 


(HI.14) 

(HI.15) 

(HI.16) 

(HI.17) 

(HI.18) 


The suppressed spacetime integrations can be restored 
by comparing these expression to the diagrams and using 
the fact that V vertices join one Higgs and two Goldstone 
lines, while Vn vertices join three Higgs lines. (These ex¬ 
pressions can be checked using the supplemental Math- 
EMATICA notebook M-) 

In terms of the SSB ansatz variables the nontrivial 
Ward identities are O (A)-scalar equations which read 


0 = Wi= / Aq^{x,z)v, (HI.19) 

V {x, y,z)v + Aq^ {x, y) - A^^ (a;, y). 

(HI.20) 


The physical meaning of these can be seen by assuming 
translation invariance, substituting 

A-^)h{p)^ f e^P-V-y^A^)^{x,y) 

Jx-y 

= p2-m^/^-SG/ff (p), (HI.21) 

and V = {—\v/3)x5V) (x — y) 5^^ (x — z)+SV where SV 
represents the loop corrections to the vertex. Matching 
powers of h (which are implicit in T,g/h and 5V) results 
in 


utog = 0, 

(HI.22) 

\v^ 2 1 r. 

— ^ + niff - Wg, = 0, 

(HI.23) 

SV (p, -p, 0) u -b Fff (p) - Eg (p) = 0, 

(HI.24) 


0 = W2 = 


















which are Goldstone’s theorem, the tree level relation 
between the particle masses, and a relation between the 
vertex correction (with one external Goldstone boson leg 
set to zero momentum) and the self-energies of the Higgs 
and Goldstone bosons respectively. The imaginary part 
of this last identity can be used to extract a relation 
between the Higgs decay rate and the off-shell Goldstone 
boson self-energy and vertex corrections. This will be 
investigated in Section IVlIl 

We now wish to impose (IHI.19I) - (IHI.20I) as constraints 
on the allowable values of f/j, A and V in the 3PIEA. 
First we review the 2PIEA case as discussed in Q, which 
imposes (IHI.19I) on P*^^) through the introduction of La¬ 
grange multiplier fields £(4 (a:), where A is an O (N) ad¬ 
joint index, and the symmetry improved effective action 
which we write in manifestly covariant form as 


It is now necessary to introduce the constraint through 
a limit process, and choose the scaling of £ in the limit 
such that the scalar equation of motion is traded for the 
constraint. To this end we set 

V f (a;,y) = (HI.32) 

Jy 

and take the limit ^ t 0. Note that, in extension of 
Pilaftsis and Teresi , one may allow separate regulators 
r]i for each Goldstone mode i = 1,...,A — 1, but there 
is nothing much to gain from this and it leads to no new 
difficulties so we take a common regulator rji = rj. m 
is an arbitrary fixed mass scale, conveniently taken to 
be ^ mH^ which serves to make 77 dimensionless. The 
modified equations of motion become 


f [<^, A, £] = r(2) [<p, ^] + y/A (^) {x) [Pt x)U • 

(HI.25) 


The transverse projector 

[PTi^,x)]^^ = Scd- 


‘Pc jx) Ifd jx) 
(x) ’ 


(HI.26) 


ensures that only the Goldstone modes are involved in the 
constraint. The equations of motion follow from Sr/S<p = 
Sf/SA = 0. 

Substituting the SSB ansatz and using translation in¬ 
variance gives 

f A,£] = r(2) [tp, A] - £>Vi, (HI.27) 


where we have absorbed group theory factors in £. The 
2PI equations of motion become 


5r(2)/i/r 

dv 

SAg (x, y) 

6Ah (x,y) 

0 




(HI.28) 

(HI.29) 

(HI.30) 

(HI.31) 


The factor of VT on the left hand side of the first equa¬ 
tion is the volume of spacetime, which we have divided 
by to give an intensive quantity. 

Now applying the constraint with v ^ 0 directly in the 
equations of motion would give zero right hand sides, re¬ 
ducing to the standard 2PI formalism. This is valid in the 
full theory because the Ward identity is satisfied. How¬ 
ever, this is impossible in the case where the 2PI effective 
action is truncated at finite loop order because the actual 
Ward identity obeyed by the 2PIEA is 
The manifestation of this fact in the symmetry improve¬ 
ment formalism is a singularity: £ —>■ 00 as v J A^^ —^ 0 
so as to leave a finite right hand side in the first equation 
of motion. 


ar(2)/yT _ 

^ 5 

OV V 

_ _£rf^ 6 
5Ag{z,w) V 


(HI.33) 

(HI.34) 


If we choose to scale 77 and the £ such that £0 = £r]/v is 
a constant and £ 77 ^/ 7 ; —?> 0 then 


ar( 2 )/i/r 

dv 

^r( 2 ) 

5Ag (z, w) 


= £ow^, 


= 0 , 


(HI.35) 

(HI.36) 


in addition to the Ward identity and the Ah equation of 
motion. In practice, in the symmetry broken phase, one 
simply discards the first equation of motion and solves 
the second one in conjunction with the Ward identity, 
which suffices to give a closed system. In the symmetric 
phase w = 0 and the Ward identity is trivial, but P^^^ 
also does not depend linearly on u, hence one can take 
the previous equations of motion with £0 = 0. Note that 
we can, and do, keep a nonzero toq in the intermediate 
stages of the computation to serve as an infrared regula¬ 
tor. 

To recap the procedure: first we define a symmetry 
improved effective action using Lagrange multipliers and 
compute the equations of motion. Second, note that the 
equations of motion are singular when the constraints 
are applied. Third, regulate the singularity by slightly 
violating the constraint. Fourth, pass to a suitable limit 
where violation of the constraint tends to zero. We re¬ 
quire the limiting procedure to be universal in the sense 
that no additional data (arbitrary forms of the Lagrange 
multiplier fields) need be introduced into the theory. 

We now extend this logic to the 3PI case. To that end 
we introduce the symmetry improved 3PIEA 


p(3) _ p(3) _|_ 


-Bf 


W, 


A(i) 

cd 


£l(x)W,^W (x) [Pt(<^,x)U 
, (111.37) 






















9 


where the second term is the same as the 2PI symmetry 
improvement term and the third term contains the ex¬ 
tended symmetry improvement. B is the new Lagrange 


multiplier and / 


W. 


Ail) 


is an arbitrary functional which 

vanishes if and only if its argument vanishes. Substitut¬ 
ing the SSB ansatz we obtain 


f (3) ^ p(3) _ ^ 


(III.38) 


The equations of motion are 




dv 


= -I- B 


Sf 


SW 2 [x, y) 


V (x,y,z), 


(III.39) 


Sf 


SAG(r,s) ^ Jxy SW 2 {x,y)' 


Aq (a;,r)Ag (s,?/), 

(111.40) 

=B f —-- A-i (a;, r) A’^ (s, y ), 

SAh (r, s) J^y SW 2 {x, y) ^ 

(111.41) 

r Sf 

SV {r, s, t) ^ 5yV2 (x, y) 


X (5^^^ {x — r) {y — s) {z — t), 


St(3) 


SVN{r,s,t) 

Wi = 0, 

W 2 = 0, 


(111.42) 

(111.43) 

(111.44) 

(111.45) 


where we already take the previous limiting procedure to 
eliminate £ and Wi. In (IIII.42I) we have inserted a factor 
of 1 = (5^^^ {z — f) for later convenience. Now we devise 

a limiting procedure such that the right hand sides of two 
of (|III.40I) . (IIII.41I) and (IIII.42I) vanish. The remaining 
equation must be chosen so that it can be replaced by the 
constraint lllll.45|) and still give a closed system. Note 
that (jIII.42ll cannot be eliminated because there is not 
enough information to reconstruct V from Aq/h using 
yV 2 - Thus we must eliminate either Aq or Ah, or else 
artificially restrict the form of V. 

We show that the desired simplification of the equa¬ 
tions of motion can be achieved without restricting V 
under the assumption that 5 f/SW 2 {x,y) is a space- 
time independent constant. Note that this is not re¬ 
quired by Poincare invariance (only the weaker condition 

6 f/SW 2 {x,y) = g {^x — yf'^ is mandated). We tem¬ 
porarily adopt this assumption without further explana¬ 
tion, though in Appendix we will show that it can be 
justified by the introduction of the d’Alembert formalism. 

Computing the left hand side of (jIII.421) using (III.191) 
and displaying only the two loop terms explicitly we ob¬ 


tain 


(5r(3) nf 

-J--- =- (N-l) 

5V{r,s,t) 2 ^ ’ 


Ixyz 


Xv 


V (x, y, z) + —6 {x-y)S{x- z) 


X Ah {x, r) Aq (y, s) Ac {z, t) + O (h^). 

(III.46) 


Without symmetry improvement one sets this quantity to 
zero, giving an equation equivalent to the one we derived 
in the previous section, (111.261) (up to a group theory 
factor, since the variables in the one case are Aab and V 
and in the other Aq, Ah and V). This equation is now 
modified by the symmetry improvement to 


-y (^ - 1) y z) + y ^ {x-y)S{x- 

X Ah {x, r) Aq (y, s) Aq (z, t) + 0 (hf) 

= [ Wiry— (y - ■5) (^ - ^) ■ 

Jxyz OW 2 [X, y) 

(III.47) 

Convolving with the inverse propagators 
gives 


-^{N-l)[v (a, b, c) + ^ 6 ia-b)Sia- c)^ +0 (h^) 



Sf 

SW 2 (x, y) 


^H (a;, a) Agi {y, b) 


Wi. 


(III.48) 


The right hand side now vanishes due to (IIII.44I) . With 
the regulator (IIII.32I) in place, we have that the sym¬ 
metry improvement term in the V equation of motion 
vanishes faster than the naive BSf / 5 W 2 scaling manifest 
in (IIII.42I) . Schematically, the right hand side scales as 
B {Sf / 6 W 2 ) mjjmGV ^ B {Sf/ 6 W 2 ) m^y'^/v. So long as 
BSf /SyV 2 does not blow up as fast as r]~‘^ as 77 —>■ 0 the 
symmetry improvement has no effect on V. 

Now we investigate the Goldstone propagator. Sub¬ 
stituting (III.191) into (IIII.40I) we find the symmetry im¬ 
proved equation of motion for A^: 


s) = ^OG (a s) - Sg {r, s) , (III.49) 

where we have defined the 3PI symmetry improved self¬ 
energy 


Eg ir,s) = 


2 i 


h{N-l) 


ST. 


+B 


Sf 


SW2 {x, y) 


SAg {r, s) 

Aq^ {x,r) Aq^ {s,y) . 

(III.50) 
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FIG. III.l. Sunset self-energy graph 


(S') 

Substituting in Fg to two loop order we find 


Eg (g s) = ^Tr [{N -f 1) Ag + Ah] (r - s) 

— ih f (v {a, b, r) + (a — r) (6 — r)^ 

Jabcd \ 3 J 

X V (c, d, s) Ah (a, c) Aq {b, d) 


-h 


2 i 


h{N -1) 

+ o (h^). 


B 


Sf 


6 W 2 (x, y) 


Aq {x,r)AQ {s,y) 

(III.51) 


The first term corresponds to the Hartree-Fock diagram 
fFigure im far left); the second term corresponds to the 
sunset diagrams (Figure IIII.ll) and the third term is the 
symmetry improvement term. The equation of motion 
for Ah can be written in the same form with a suitable 
definition for a symmetry improved self energy T,h (r, s), 
where the symmetry improvement term now has the form 
^Bj{Sf/6W2)A-H^A-^\ 

If we assume Sf /SyV 2 is constant, we find that Eg and 
Eh scale as {B 5 f/ 5 W 2 ) rriQ ^ {B 5 f/ 5 W 2 ) 

{B5f /SyV 2 ) rn^H ^ {BSf/ 5 yV 2 )Tri\jrf’ respectively. Thus, 
by choosing a regulator such that BSf / 6 yV 2 goes to a fi¬ 
nite limit, the equations of motion for V and Ag are un¬ 
modified and the equation of motion for Ah is modified 
by a finite term. This is the desired limiting procedure. 
Adopting it gives the final set of equations of motion: 


5AG{r,s) 

(HI.52) 

SV{r,s,t) 

(HI.53) 

«r«) ^ 

SVn {r, s, t) 

(HI.54) 

Wi = 0, 

(HI.55) 

W 2 = 0. 

(HI.56) 


IV. RENORMALIZATION 

Here we undertake a general description of the renor¬ 
malization problem at zero temperature. Our detailed 
considerations follow in Sections IIV Al and IIVBI for two 
and three loop truncations respectively. Finite tempera¬ 
ture results are given for the Hartree-Fock approximation 
in Section |Vl The two loop renormalization of the theory 
in Section IIV Al is non-trivial already even though the 
vertex equation of motion can be solved trivially. This is 


because the symmetry improvement breaks the nPIEA 
equivalence hierarchy by modifying the Higgs equation 
of motion. 

Generically, modifications of the equations of motion 
following from the 2PIEA will lead to an inconsistency 
of the renormalization procedure since the 2PIEA is self- 
consistently complete at two loop order (in the action, 
i.e. one loop order in the equations of motion). However, 
we will see that the wavefunction and propagator renor¬ 
malization constants (normally trivial in at one loop) 
provide the extra freedom required to obtain consistency. 
Then in Section IIVBI we will renormalize the theory at 
three loops. Non-perturbative counter-term calculations 
are generally much more difficult than the analogous per¬ 
turbative calculations, hence many of the manipulations 
were performed in a supplemental Mathematica note¬ 
book [ 1 ^. The results of this section are finite equa¬ 
tions of motion for renormalized quantities which must 
be solved numerically. We leave the numerical implemen¬ 
tation to future work, except in the case of the Hartree- 
Fock approximation. 

We wish to demonstrate the renormalizability of the 
equations of motion (IHI.52|) - (IHI.56I) . First we examine 
the symmetric phase, since on physical grounds SSB is 
irrelevant to renormalizability. In the symmetric phase 
?; = 0 and the Ward identity (IHI.55I) is trivially satisfied, 
while (IHI.56I) requires Ag = Ah as expected. Further, 
iteration shows that (IHI.53I) and (IHI.54I) have the solu¬ 
tion V = V/v = 0 as expected on general grounds: there 
is no three point vertex in the symmetric phase. As a 
result, the symmetry improved 3PIEA in the symmetric 
phase is equivalent to the ordinary 2PIEA 

f (3) = 0 , A, V = 0] = [if = 0, A], (IV.l) 

which is known to be renormalizable, either by an implicit 
construction involving Bethe-Salpeter integral equations 
or an explicit algebraic BPHZ (Bogoliubov-Parasiuk- 
Hepp-Zimmerman) style construction which has non¬ 
trivial consistency requirements (but which has been 
shown to be equivalent to the Bethe-Salpeter method) 
[IS US mi- Thus, only divergences arising from nonzero 
V pose any new conceptual problems. 

We will extend the BPHZ style procedure of [10, 
which was adapted to symmetry improved 2PIEA by 
[S and to 3PIEA for three dimensionalpure glue QCD 
(without symmetry improvement) by (Sj. The essence 
of the procedure is quite simple. Consider for example 
the quadratically divergent integral f iAQ/nid)- Since 
Aq/h{<i) is determined self-consistently this is a com¬ 
plicated integral which must be evaluated numerically. 
However, the UV behavior of the propagator should ap¬ 
proach q~'^ as g —>■ 00 . (Note that Weinberg’s theorem 
[10] implies that the self-consistent propagators have this 
form up to powers of logarithms [22| , though renormal¬ 
ization group theory shows the true large-momentum be¬ 
havior of the propagators is a power law with an anoma¬ 
lous dimension. This implies that a truncated nPIEA 
does not effect a resummation of large logarithms.) Now 
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we can add and subtract an integral with the same UV 
asymptotics: 


i^G/H (q) 



i 

— fj? -\- ie 


(IV.2) 


where is an arbitrary mass subtraction scale (not a 
cutoff scale). The first term is now only logarithmically 
divergent and the second term can be evaluated analyt¬ 
ically in a chosen regularization scheme such as dimen¬ 
sional regularization. A further subtraction of this kind 
can render the first term finite. 

We write the renormalized propagators as 

~ ~ ^G/H (p ), (IV.3) 

^G/H ip) = ^G/H iP) + ^%/H iP) + ^h/H iP) , (IV.4) 


where me /h is the physical mass and the (renormalized) 
self-energies have been separated into pieces according 
to their asymptotic behavior: (p) ^ (lnp)‘^\ 


yO 
^G/H 


(Inp)'"" and 


' p ^ as p 


oo respectively. 


The pole condition requires ^q/h {p^ = ~ 

We also introduce the auxiliary propagator = 

( 


P P ^G/H 
expanded in 


. The propagator Aq/h can be 


^g/h{p) = ^G/niP) 


+ 

+ 0 


A^ 


^G/h(p)] (’ 


m; 


''G/H 


- P 


■ '^G/H 


' ^g/h) 




1 3 


This allows us to extract the leading order asymptotics 
of diagrams as p —>■ oo. 

We now do a similar analysis to isolate the lead¬ 
ing asymptotics for V at large momentum. Sup¬ 
pressing O (N) indices we can write V (pi,p 2 ,P 3 ) = 
(v’ where pi -f P 2 + Pa = 0. Now V ^ 0 

as V ^0 implies that f {xi,X2,X3) ^ x“ (lnx)''^ where 
a < 1 and y is representative of the largest scale among 
Xi, X 2 and X 3 . Now consider the vertex equation of mo¬ 
tion (III.26P or Figure III. 21 The triangle graph goes like 

/^^3a-6(lj^£)3c3-3ci _ ^3a-2 y^)3c3-3ci ^ ^ 2/3 

or (Inx)^'*'^'^^”^'^^ if a = 2/3, which is dominated by 
the bubble graph which goes like ^ 

X“ (Inx)'^^”^'^^ if a ^ 0 or (Inx)^'*'^^”^'^^ if a = 0. Thus 
to a leading approximation the large momentum behav¬ 
ior is obtained by dropping the triangle graph from the 
equation of motion. This can also be seen by taking 
u —>■ 0 at fixed pi which suppresses the triangle graph 
relative to the bubble graph. 

We now define auxiliary vertex functions and 
which have the same asymptotic behavior as V and Vn 



/ \ / \ / \ 


FIG. IV.1. Defining equation for the auxiliary vertex func¬ 
tions V'^ and V/), obtained by taking the corresponding equa¬ 
tions of motion for V and Vjv, dropping the triangle graphs, 
and making the replacements V —^ Viv —>■ V^, and 

Aq/h Aq^u- The filled triangle represents the auxil¬ 
iary vertices and the lines represent the auxiliary propagators. 
Crossed circles are bare vertices as before. 


respectively, though depend only on the auxiliary prop¬ 
agators. We define and by taking the equations 
of motion for V and Fv, dropping the triangle graphs, 
and making the replacements F —>■ F'^, Fat —>■ Fjy, and 
Ag/h Aqjij. These equations are shown in Figure 
IIV.II This gives a pair of coupled linear integral equa¬ 
tions, analogous to the Bethe-Salpeter equations, for F^ 
and which can be solved explicitly by iteration. (De¬ 
tails of this calculation are presented in Section lIVBp . 
Unfortunately the result is only analytically tractable in 
fewer than 1-1-3 dimensions, so we confine the analytical 
results depending on the explicit forms of and F/( to 
this case. In the physically most interesting case of 1-1-3 
dimensions, F^ and F/^ must be numerically determined 
at the same time as F and Fat. 

By using these auxiliary propagators and vertices we 
can isolate the divergent contributions to the equations 
of motion and so obtain the required set of counter-terms 
to remove them. 


A. Two loop truncation 

The theory simplifies dramatically at two loop order. 
It follows from (III.261) that at this order F —>■ Vb (up 
to a renormalization). Substituting this into the action 
gives, apart from the symmetry improvement terms, the 
standard 2PIEA. This is an example of the equivalence 
hierarchy previously discussed. Another simplification 
is that the logarithmic enhancement of the propagators 
in the UV due to '£‘g/h vanishes at this level C^g/h I® 
generated by the diagram < 1)5 appearing at three loop 
order). In this case Aq = = A^ = (jp — 

However, the reduction is not trivial because now the 
Higgs equation of motion has been replaced by a Ward 
identity. The equations of motion reduce to 


Ag^ {x,y) 



O'" + m'^ + 


{x - y) 
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if) 

“ + 1) {x, x) {x - y) 

o 
if) 

- — AAh {x, x) (x - y) 

6 

-^AVAj^(x,y)AG(x,y), (IV.6) 

\,,2 

(x, y) = — {x-y) + Ag^ (x, y), (IV.7) 

VTO^ = 0. (IV.8) 

The first line is the tree level term, the second and third 
lines are the Hartree-Fock self-energies, the fourth line 
is the sunset self-energy, and the last two lines are the 
Ward identities VV 2 and Wi respectively. 

To renormalize the theory we regard all parameters 
heretofore as bare parameters and introduce renormal¬ 
ized counterparts using the same letters: 


',v) -)> , 

(IV.9) 

—>• Z~^Z^ (rn? 8m^) , 

(IV.IO) 

A^ Z-2 {X + 8X), 

(IV.ll) 

A -)> ZZaA, 

(IV.12) 

V Z-^l'^ZvY. 

(IV.13) 


Hereafter whenever we refer to a bare parameter we in¬ 
dicate this with a subscript “B,” e.g. etc. The wave- 
function renormalizations for A and V can be obtained 
from their definitions A ^ and AAAH ^ {<P4‘4>) 
respectively. Due to the presence of composite opera¬ 
tors in the effective action, additional counter-terms are 
required compared to the standard perturbation theory: 
Sttiq and dAo for terms in the bare action, dml for one 
loop terms, dA^ for terms of the form (picjiiAjj, dXf for 
(jiicjijAij terms, 8 X 2 for AuAjj and 8 X 2 for AijAij. Sim¬ 
ilarly, A and V are given independent renormalization 
constants Za and Zy respectively. We give terms in the 
sunset graphs a universal 8 X counter-term. 

The renormalized equations of motion are (see Ap¬ 
pendix |B] for more detail) 


(p) — — w? — 8m\ — Za 

6 

- ^ [(A + 1) A + (A - 1) 8 X^ + 2dAf ] ZITg 


- - (A + 8 X^) ZITh 


ih 


(A -I- 8 X) V 


Z%Ihg {p) ■. 


A-i r ^ 7 + , A-i / ^ 

{p) = -Za -^- + Aq {p), 

vrriQ = 0, 


(IV.14) 

(IV.15) 

(IV.16) 


where for convenience we have defined 


Tg/h = 


i^G/H (q) ) 


(IV. 17) 


Zhg {p) = / (<?) iAg (p-q)- (IV. 18) 

Jq 

The tadpole integrals Tg/h correspond to the Hartree- 
Fock graphs and Ihg (p) corresponds to the sunset self¬ 
energy graph. From these we identify the self-energy 
parts 


{p) = {ZZa-1)p^ = 0, 


-‘G/H (p) + Z'g/H (p) — 


(A -|- i5A) V 


(IV.19) 

Z^Ihg ip) ■ 
(IV.20) 


Note that the Goldstone and Higgs self-energies are equal 
to this order as a consequence of the vertex Ward identity. 
This is essentially where our treatment differs from Q . 

By using the auxiliary propagators to extract the di¬ 
vergences in Tg/h and Ihg (p) and absorbing them into 
the counter-terms (see Appendix [g we find the finite 
equations of motion 


(p) =p^-m^-^v^-^{N + 1) A7^" - ^AT^" 

+ [zfw (p) - Z^G (mG)] , (IV.21) 

(IV.22) 

The finite parts and if^Q (p) are 


I%ip)=lHG{p)-I^. (IV.23) 

= 'Tg/h -r^ + i [rnl/H - 

- /*[A'^(g)]"E'^(g), (IV.24) 

Jq 

where the auxiliary quantities are 

= /* (g), (IV.25) 

J q 

I^= f [iA>^ (g)]" , (IV.26) 

Jq 


iA>^ {l)iA^ {q + ^)-I^^ . 

(IV.27) 

(For details see Appendix [BJ) These equations are the 
main result of this section. We expect they could 
be solved numerically using techniques similar to 
though we leave the numerical implementation for later 
work. 


S'" ( 9 ) = -ih j 
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B. Three loop truncation 

We consider now the three loop truncation of the ef¬ 
fective action. The vertex equation of motion is shown 
in Figure HL^ and we have already argued that the lead¬ 
ing asymptotics at large momentum are captured by the 
auxiliary vertex defined by its equation of motion in Fig¬ 
ure IIV.II Subtracting these two equations we find that 
the right hand side is finite (indeed the auxiliary vertices 
were constructed to guarantee this). Thus the problem 
of renormalizing the vertex equation of motion reduces to 
the problem of renormalizing the auxiliary vertex equa¬ 
tion of motion. 

It is temporarily more convenient to go back to the 
O {N) covariant form we had before introducing the SSB 
ansatz. Introduce the covariant auxiliary vertex 
which is related to and by an equation analo¬ 
gous to (IIII.I3I) . Iterating the equation of motion we find 
the solution 


^abc ~ ^abcdefVQdef ^ (IV.28) 


where the six point kernel Kabcdef obeys the Bethe- 
Salpeter like equation 


^abcdef — ^ad^be^cf 


3! 


E - 


3ih 


^'K(a)h 


X (IV.29) 


where is a sum over permutations of the incoming 
legs. This equation is shown in Figure ITV.21 (I IV. 291) can 
be written in a form that makes explicit all divergences 
(see Appendix [C]) and replaces the bare vertex IV by a 
four point kernel ~ A/ (1 -I- AI^). 

In fewer than four dimensions is finite and every 

correction to the tree level value is asymptotically sub¬ 
dominant. Thus the leading term at large momentum is 
the tree level term and, instead of the full auxiliary vertex 
as we have defined it, one can simply take = Voa&c, 
dramatically simplifying the renormalization theory. A 
similar simplification happens to the auxiliary propaga¬ 
tor due to the logarithmic (rather than quadratic) diver¬ 
gence of ^s-generated self-energy in < I -|- 3 dimensions. 
This confirms statements made in the literature (sup¬ 
ported by numerical evidence though without proof, to 
our knowledge) to the effect that the asymptotic behavior 
of Green functions is free (e.g. Q). 

I 


Unfortunately, the situation is much more difficult in 
four dimensions and the renormalization of the nPIEA 
for n > 3 in d > 3 remains an open problem, both in 
general and in the present case. The problem can be 
seen from the behaviour of the auxiliary vertex which is 
discussed further in Appendix [C] For the sake of obtain¬ 
ing analytical results we restrict the rest of this section 
to < I -f 3 dimensions. The renormalization of the I -I- 3 
dimensional case is left to future work. 

We derive the counter-terms for I-I-2 dimensions in Ap¬ 
pendix |D] The are only two interesting comments about 
this derivation: the first is that we require an additional 
(non-universal) counter-term for the sunset graph lin¬ 
ear in V] the second is that, consistent with the super- 
renormalizability of 0'^ theory in I -|- 2 dimensions, only 
(5mf is required to UV-renormalize the theory. Every 
other counter-term is finite and exists solely to maintain 
the pole condition for the Higgs propagator despite the 
vertex Ward identity. The resulting finite equations of 
motion are 


= - ( dud^ -I- -I- 


A 2 
6^ 

- [^G (P) - ( wg )] , 

for the Goldstone propagator. 


(IV.30) 


V = - 


Xv 


3 

ih 

ihX r 
“6 


Vn (V) (A//)^ Ag-I-(U) AniAcY 


Vn {Ah)^ -I- (V -I-1) U (Ag)^ -I- AV Ac Ah 


for the Higgs-Goldstone-Goldstone vertex, and 


(IV.31) 


Vat = —Xv 


ih 

ihX 


{N - 1) {Vf {Acf + {VNf [AhY 


{N-l)VAGAG + WN{AHf , (IV.32) 


for the triple Higgs vertex. 

The finite Goldstone self-energy is 


(p) = - (TV + 1) A (Tg - rn - ^A {Th - rn 


6 


— ih 




6 

AhAgV + h^ 


Vn {V)^ (Ah)^ (Ag)'^ + (V) AhAh{AgY 


— VVn {Anf Ag + {N + I) VVAh (Acf + 3UU (Ag)" A^ A^^ 



























14 


+ —rTT 


{N + l){AGf + AhAhAg 



FIG. IV.2. Solution for the auxiliary vertex function in terms 
of a six point kernel ICabcdef which is represented by the blue 
box (the indices run from top to bottom down the left side, 
then the right). The vertical black bar in the kernel equa¬ 
tion of motion represents symmetrization of the external lines 
(with a factor of 1/3!). 




' I \ 


A N 

/ \ ' 
- * - - » 



FIG. IV.3. Feynman graph topologies appearing in the self¬ 
energy function Eg (p) in (IIV.33II . 


where the bball integral is = 

A^ {q) A^^ (p) A^^ {p + q). The graph topologies 
are shown in Figure HV.31 Finally, the Higgs equation of 
motion is 


(P) = (™G + Sg - mjj) 


X 


_ V (p, -P, 0) 

V (niff, 


+ ^G^ (P) ■ 


(IV.34) 


The unusual form of this equation is a result of the pole 
condition A]j^ (mH) = 0. We defer the numerical imple¬ 
mentation of these equations to future work. 


V. SOLUTION OF THE HARTREE-FOCK 
APPROXIMATION 

In the Hartree-Fock approximation one drops the 
Ihg (p) term in the two loop equations of motion, or 
equivalently drops the sunset diagram. In this case the 
problem simplifies dramatically because the self-energy 
is momentum independent. The machinery of the aux¬ 
iliary propagators introduced previously is now unneces¬ 
sary and Tg/h = H written as the sum 


-{N + 2)B>^ , 


(IV.33) 


of divergent and finite parts which can be evaluated in 
closed form. In the Matsubara formalism at finite tem¬ 
perature T the time contour is taken on the imaginary 
axis with periodic boundary conditions of period —i^-, 
where /3 = 1/T. Integration over the timelike momen¬ 
tum component p^ becomes a sum over discrete Matsub¬ 
ara frequencies = 27rn//3, n = 0, ±1,±2,---. Using 
standard tricks [l^ the sum over frequencies can be per¬ 
formed, giving 


'T-'OO 

'G/H 


'G/H 


^7—vac 
'G/H 


^ 7 -th 

'G/H 


"‘'G/H 

1071^ 


— 7 -f 1 -b In ( 47 r) 




'G/H + 'G/ffi 


'"‘G/H 

f 1 


In 


'-G/H 


/q Wq - 1 ’ 


(V.l) 

(V.2) 

(V.3) 

(V.4) 


where the divergent and finite vacuum parts have been 
evaluated using MS in d = 4 — 2e dimensions at the 
renormalization point p (note that Q adopt a slightly 
different convention for TqJ/j which amounts to a redef¬ 
inition of p not affecting physical results). 7 ~ 0.577 is 
the Euler-Mascheroni constant. In the thermal part q 
is the spatial momentum vector and ujq = ^-b 

We substitute these expressions into the equations of mo¬ 
tion (IIV.14p - (IIV.16|) and demand that the kinematically 
distinct divergences proportional to v, 7^/// inde¬ 

pendently vanish. The other renormalization conditions 
are that residue of the pole of the propagator equals one, 
which requires ZZ^ = 1 , and that the tree level rela¬ 
tion m'^ = ^ + rriQ holds at zero temperature. These 
conditions determine the renormalization constants 


Z = Za = 1, 


2 + 2 ) hXrn? (eK + 1 ) 

om^ = -- 


9071^ 


7 _ fiA(jV-|-2)(eK+l) ’ 
967r^€ 


cxA c 2 

[N + 2) 

SX^ = SX^ = ^^6Xf, 
^ ^ N + A 


(V.5) 

(V.5) 

(V.7) 

(V. 8 ) 


where k = 1 — 7 -bln 47 r« 2.95. Note that the undeter¬ 
mined constant dA can be consistently set to zero at this 
order. The finite equations of motion are 


ml = ^ (W + 1) 73” + ^T^”. 

2 2 
ruH — b rriQ, 

vtUq = 0. 


(V.9) 

(V.IO) 

(V.ll) 
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Finally, if we demand the zero temperature tree level 
relation {T — 0) = v'^ = —6m^/A we must set the 
renormalization point /i^ = ffij^ = mjj (T = 0) = An^/3. 

The analogue of the equations of motion (jV.9|) - (jV.lll) 
corresponding to previous work on the symmetry im¬ 
proved 2PIEA is ([^, [13| generalized to arbitrary N) 

= m2 + ^ (TV + 1) TJ" + (V.12) 

6 6 6 

= + + + (V.13) 

vrriQ = 0. (V.14) 



T (MeV) 


Note that only the Higgs equation of motion differs, as 
expected. In the standard formalism without symmetry 
improvement one replaces (IV.141) with 

0 = v(^m^ + ^v^ + f(N- l)rj" + ■ 

(V.15) 

These equations of motion, or gap equations, possess 
a phase transition and a critical point where = 
= 0. Using the result for massless particles 
'Tq/h {t^g/h = O) = T^/12, we find the same value of 
the critical temperature 


independent of the formalism used. However, the order of 
the phase transition differs in the three cases. This stands 
in contrast to the large-N approximation, which correctly 
determines the order of the phase transition but gives a 
critical temperature larger by a factor of \/3j2-\-0 {N~^) 
(see Hi,if). 

We present numerical solutions of equations (IV.911 - 
(IV.15|) with TV = 4, n = 93 MeV and friH = 500 MeV. 
These values are chosen to represent the low energy 
mesonic sector of QCD, and to enable direct compari¬ 
son with Hi]. Our results are also closely comparable 
with though they take ffiH ~ 600 MeV. The solu¬ 
tion is implemented in Python as an iterative root finder 
based on scipy. optimize.root j2^ with an estimated Jaco¬ 
bian or, if that fails to converge, a direct iteration of the 
gap equations. The Bose-Einstein integrals in (IV. 4p can 
be precomputed to save time. We show the results for 
the scalar field v, Higgs mass uih and Goldstone mass 
me in Figures rvTTl IV.21 and IV. 31 respectively. 

Figure [VT] shows v{T), the order parameter of the 
phase transition. Below the critical temperature there 
is a broken phase with v 0, but the symmetry is 
restored when n = 0 above the critical temperature. 
Note, however, that the unimproved and symmetry im¬ 
proved 3PIEA have unphysical metastable broken phases 
at T > T*, signalling a first order phase transition. 
The symmetry improved 2PIEA correctly predicts the 
second order nature of the phase transition. Though 


FIG. V.l. Expectation value of the scalar field v = {(j)) as 
a function of temperature T computed in the Hartree-Fock 
approximation using the unimproved 2PIEA (solid black), the 
Piiaftsis and Teresi symmetry improved 2PIEA (dash dotted 
blue) and our symmetry improved 3PIEA (solid green). In 
the symmetric phase (dashed black) all methods agree. The 
vertical grey line at T « 131.5 MeV corresponds to the critical 
temperature which is the same in all methods. 


unphysical, the symmetry improved 3PIEA behavior is 
much more reasonable than the unimproved 2PIEA: the 
strength of the first order phase transition is reduced 
and the metastable phase ceases to exist at a tempera¬ 
ture much closer to the critical temperature than for the 
unimproved 2PIEA. Figure IV.21 shows the Higgs mass 
mn (T). The phase transition behavior above is seen 
again, and again all three methods agree in the symmet¬ 
ric phase, giving the usual thermal mass effect. Finally, 
Figure IVTSl shows the Goldstone boson mass. The unim¬ 
proved 2PIEA strongly violates the Goldstone theorem, 
but both symmetry improvement methods satisfy it as 
expected. Note that the Goldstone theorem is even sat¬ 
isfied in the unphysical metastable phase predicted by the 
symmetry improved 3PIEA. All three methods correctly 
predict tog = rnn in the symmetric phase. 


VI. TWO DIMENSIONS AND THE 
COLEMAN-MERMIN-WAGNER THEOREM 


Recall that the Coleman-Mermin-Wagner theorem 
27j |. which has been interpreted as a breakdown of the 
Goldstone theorem [28|, is a general result stating that 
the spontaneous breaking of a continuous symmetry is 
impossible ind = 2ord=l-|-l dimensions. This oc¬ 
curs due to the infrared divergence of the massless scalar 
propagator in two dimensions. We show that the symme¬ 
try improved gap equations satisfy this theorem despite 
the direct imposition of Goldstone’s theorem. Thus sym¬ 
metry improvement passes another test that any robust 
quantum field theoretical method must satisfy. (Note 
that symmetry improvement is not required to obtain 
consistency of nPIEA with the Coleman-Mermin-Wagner 
theorem, but neither does it ruin it.) 
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FIG. V.2. The Higgs mass rriH as a function of tempera- 
tnre T computed in the Hartree-Fock approximation nsing 
the nnimproved 2PIEA (solid black), the Pilaftsis and Teresi 
symmetry improved 2PIEA (dash dotted blue) and our sym¬ 
metry improved 3PIEA (solid green). In the symmetric phase 
(dashed black) all methods agree. The vertical grey line at 
T « 131.5 MeV corresponds to the critical temperatnre which 
is the same in all methods. 


equations possess no solution for rriQ = 0. It is clear that 
if u 0, diverges as ttiq —>• 0 if we take fi as a con¬ 
stant, and diverges if we take oc uiq as tuq —>■ 0 . 
Either way there is no solution. At finite temperature 
the Bose-Einstein integral also has an infrared di¬ 
vergence as TTia 0 which does not cancel against the 
singularity of the vacuum term. It can be shown that the 
singularity is due to the Matsubara zero mode. 

For u = 0 on the other hand, the gap equations reduce 
to 

rn% = ml. = m^ - + 2)X\n(^^'^ , (VI.2) 

which always has a positive solution. If > 0 then one 
can choose the renormalization point = uiq so that 
the tree level relationship = rri^ holds. If < 0 a 
positive mass is dynamically generated and one requires 

... .9 9 / 24-771 I \ 

a renormalization point fj. > rriQ exp I ) ^lon- 

perturbatively large in the ratio A/ |m^|, reflecting the 
fact that perturbation theory is bound to fail in this case. 



VII. OPTICAL THEOREM AND DISPERSION 
RELATIONS 

In this section we examine the analytic structure of 
propagators and self-energies in the symmetry improved 
3PI formalism. A physical quantity of particular interest 
is the decay width E// of the Higgs, which is dominated 
by decays to two Goldstones. Fh is given by the optical 
theorem in terms of the ima gina ry part of the self-energy 
evaluated on-shell (see, e.g. Chapter 7): 

- mBFH = IrnSn (mn) ■ (VII.1) 


FIG. V.3. The Goldstone mass ma as a function of temper¬ 
ature T computed in the Hartree-Fock approximation using 
the unimproved 2PIEA (solid black), the Pilaftsis and Teresi 
symmetry improved 2PIEA (dash dotted blue) and our sym¬ 
metry improved 3PIEA (solid green). In the symmetric phase 
(dashed black) all methods agree. The vertical grey line at 
T « 131.5 MeV corresponds to the critical temperature which 
is the same in all methods. 


The general statement of the result is that , S (x, 0 ) 
diverges whenever massless particles appear in loops in 
d = 2, thus, by (III.171) and (IIII.19I) u = 0 and a mass 
gap is generated. We will show this explicitly using the 
Hartree-Fock gap equations (IV.9l) - (IV.llll . where in two 
dimensions 


•C'“(MS) = -Tl„(^^y (VI.l) 

(Note that the renormalization can be carried through 
without difficulty in two dimensions. Only the Smf 
counter-term is needed.) We must show that the gap 


(This is valid so long as Tjj <C tub , otherwise the full en¬ 
ergy dependence of Sb (p) must be taken into account.) 
The standard one loop perturbative result gives 


V-1 n 

2 IGniTiB \ 3 y ’ 


(VH.2) 


which comes entirely from the Goldstone loop sunset 
graph. Each part of this expression has a simple inter¬ 
pretation in relation to the tree level decay graph (Figure 
IVH.ll) . The V — 1 is due to the sum over final state Gold¬ 
stone flavours, the factor of 1/2 is due to the Bose statis¬ 
tics of the two particles in the final state, the h/WirmB 
is due to the final state phase space integration and the 
(Au/3)^ is the absolute square of the invariant decay am¬ 
plitude. 

The Hartree-Fock approximation fails to reproduce 
this result regardless of the use or not of symmetry 
improvement. This is because there is no self-energy 
apart from a mass correction. Thus the Hartree-Fock 
approximation always predicts that the Higgs is sta¬ 
ble. Attempts to repair the Hartree-Fock approximation 
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through the use of an external propagator lead to a non¬ 
zero but still incorrect result. This is because an un¬ 
physical value of me still appears in loops. Satisfactory 
results are obtained within the symmetry iimroved 2PI 
formalism for both on- and off-shell Higgs Q. Here we 
show that the symmetry improved 3PIEA can not yield 
a satisfactory value for P// at the two loop level. 

From (jIV.22(l 


ImS// (p) = ^{N +1) AlmTc" + 

-^(y) Im[*l|G(p)] 


= ^ (TV + 1 ) AlmTc + ^AlmTff 
6 6 
/ \ 2 


(VH.3) 


which can be written in terms of the un-subtracted Tq/h 
and Ihg because all of the subtractions are manifestly 
real. Now we show that ImTc/ff = 0. To do this we 
introduce the Kallen-Lehmann spectral representation of 
the propagators [s^ 

Ac/ff (g) = r ds f , (VH.4) 

Jo q — s le 

where the spectral densities pc /h (s) are real and positive 
for s > 0 and obey the sum rule 


where the last equality holds at two loop order (we have 
adapted the standard formula to our renormalization 
scheme). 

(Note that this standard theory actually conflicts with 
the asymptotic (inp^) ^ form assumed for the self¬ 
energy when the 4)5 graph is included, so that our ar¬ 
gument must be refined at the three loop level. The es¬ 
sential problem is that the self-consistent nPI propagator 
is not resumming large logarithms. However, it seems un¬ 
likely that a refinement of the argument to account for 
this fact will change the qualitative conclusions of this 
section since, as will be shown shortly, the predicted Pj/ 
is wrong by group theory factors in addition to the 0 ( 1 ) 
factors which could be compensated by a modification of 

PG/H-) 

Then 

Im f i f = Im / dp^PG/H (p^) = 0 . 

Jq Jo q — p + le Jq 

(VH. 6 ) 

This allows us to obtain a dispersion relation relating the 
real and imaginary parts of the self-energies 


0 = Im 


-’rnc/H -^G/h {q) 


2 2 
q — m 


G/H ~ RsSg/jt (g) 


nT^G/H 


(g) -b [imEc/ff (g)] 

(VH.7) 


dspG/H (s) = (ZZa) ^ = 1, 


Finally we have left to compute Im [Hhg (p)] which 
(VH.5) can be written 

I 


ipN (si) 


ipG (S 2 ) 


p poo poo 

Im [iJjfG (p)] = Imz / / dsi / ds 2 -„ , . , ^ 

Jq Jo Jo g^ — Si + {p — q) — S 2 -b ie 


= lmi / dsi / ds 2 / 9 Ar(si)pG(s 2 ) 

Jo Jo Jq 

1 


q^ — Si + ie [p — g)^ _ -b ie 


IGtt^ 


pOO pOO pi / 

/ dsi / ds 2 Piv (si)pG (S 2 )lm / dxlni 

Jo Jo Jo \ 


—X (I — a:) -b xsi -b (I — x) S 2 — 


. (VH. 8 ) 


The imaginary part of the x integral is only nonzero 
for yfsi -b < \/qZ. We denote the region of si ,2 inte¬ 
gration by H. Then the imaginary part of the x integral 
can be evaluated straightforwardly, giving 


Im [Hhg (p)] = 


dsids2Piv (si)pG (S 2 ) 


167rp2 jsj 

X \Jp^ - (v/sT -b \/si')^ 


(VH.9) 


Now, since the each term of the integrand is positive 
and the square root is < p^ we have 

If I 

Im [iThg (p)] < / dsids 2 Piv (si) Pg (S 2 ) < 7 ^, 

IbTT IbTT 

(VH.IO) 

using the sum rule for pm/g (s). Finally we have 


Th < 


h 

IbTTTOjf 



(VH.ll) 


which is smaller than the expected value for all > 3. 
N = 2 and 3 are cases where one could possibly obtain 
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an (accidentally) reasonable result, depending on the pre¬ 
cise form of the spectral functions, but it is clear that one 
should not generically expect a correct prediction of F// 
from the symmetry improved 3PIEA at two loops. The 
source of the problem is the derivation of the two loop 
truncation where we dropped the vertex correction term 
in (|IIL24I) . resulting in a truncation of the true Ward 
identity (|IIL20p that keeps the one loop graphs in Eg 
but not in V. The diagram contributing to F^f above 
is thus the Goldstone self-energy shown in Figure IVn.2l 
which has the incorrect kinematics and lacks the required 
group theory (N — 1) and Bose symmetry (1/2) factors 
as well. In fact, a perturbative evaluation of Figure lVII.2l 
gives Fjj = 0 due to the threshold at What 

we have shown is that no matter the form of the exact 
spectral functions, there cannot be a non-perturbative 
enhancement of this graph large enough to give the cor¬ 
rect F// for > 3. The neglected vertex corrections 
give a leading O {K) contribution to F^f which must be 
included. 

Now we consider the three loop truncation of the sym¬ 
metry improved 3PIFA. Since one loop vertex corrections 
appear at this order we expect that F h should be correct 
at least to O (h). Since the previous result was incorrect 
by group theory factors already at O (H) our task sim¬ 
plifies to seeking only the O (h) decay width, and so we 
make use of only the one loop terms in the Higgs equa¬ 
tion of motion, which are displayed in Figure lVII.31 Fur¬ 
thermore, by iterating the equations of motion we may 
replace all propagators and vertices by their perturba¬ 
tive values to O (h). This will leave out contributions of 
higher order decay processes such as iF —>■ GGGG. We 
leave the numerical task of computing the exact decay 
width predicted by the symmetry improved 3PIEA to 
future work. 

The contributions of the various terms in Figure IVII.31 
to the imaginary part of Ejj can be determined using 
Cutkosky cutting rules [1^. In particular, the Hartree- 
Fock diagram and the first bubble vertex correction di¬ 
agram (left diagram, bottom row Figure IVII.31) have no 
cuts such that all cut lines can be put on shell. Also, 
cuts through intermediate states with both Goldstone 
and Higgs lines contribute to the process H ^ HG, 
which vanishes due to the zero phase space at thresh¬ 
old. This means we can drop the sunset diagram and the 
last bubble vertex correction (right diagram, bottom row 
Figure rVH.3l) . Similarly cuts through two intermediate 
Higgs lines can be dropped since H ^ HH is impossi¬ 
ble on shell. This mean we can drop the contributions 
to the triangle and remaining bubble diagram where the 
leftmost vertex is Vn rather than V. The contributions 
we are interested in can now be displayed explicitly: 


-T,h D Vv 


D V 


ih 



1 

{£ — pf' — uiq + ie 


I 1 

P — toq — rnjj -I- ie 

ih\ f \v\ 

+ -j-(« + i)(-Yj 

^ L{£-pf -m% + ie£^ -rri^G + ie 

where the first and second term are the triangle and bub¬ 
ble graphs respectively. We now cut the Goldstone lines 
by replacing each cut propagator — rriQ + ie) 

—2niS (p^ — ttiq) to give —2ilmT,H (because the cutting 
rules give the discontinuity of the diagram, which is 2i 
times the imaginary part), yielding 


—2HmS^f D —ihv 


Xv 


A 


3 / ® 


+ -(7V + 1) - 


Xv 


X J 2ttS ^(£ — p)^J 27rS (£^) 

^ “ 1 ) J ((^ - pf) 2705 , 

(VH.13) 


The integral can be evaluated by elementary tech¬ 
niques, giving 


J^2ttS (^{£ - pf) 2TTd (£2) = ^ 

X 6 {£^ 




d^F 

2£-p+p'^) d (£^) 
d£odl4:'Kp 


X 5 {-2£omH + rrfj) 5 (Fq - 


_ 1 
Stt ’ 


I 5(^-0 


2 mH 2^ 


(VH.14) 


and finally 


= + 0 ( 6 »). 

(VH.15) 

This exactly matches the expected Fjj, including group 
theory and Bose symmetry factors. The full non- 
perturbative solution will give corrections to this ac¬ 
counting for loop corrections as well as cascade decay 
processes H GG —>■ {GGf —>■•••. We leave the eval¬ 
uation of this to future work, however, we have shown 
that the one loop vertex corrections are required to get 
the correct F// at leading order. 
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FIG. VII.l. Tree level of decay of the Higgs {H) to two Gold- 
stone bosons a = 1, • • • , — 1. 


a 



H 


FIG. VII.2. The self-energy diagram from (IIV.22II which, due 
to the inconsistent truncation of the Ward identity, gives the 
incorrect absorptive part to the Higgs propagator in the two 
loop truncated symmetry improved 3PIEA. a = 1, • • • ,N — 1 
labels Goldstone boson lines and H labels the Higgs boson 
line. 


VIII. DISCUSSION 

The symmetry improvement formalism of Pilaftsis and 
Teresi is able to enforce the preservation of global sym¬ 
metries in two particle irreducible effective actions, allow¬ 
ing among other things the accurate description of phase 
transitions in strongly coupled theories using numerical 
methods that are relatively cheap compared to lattice 
methods. As an example of this, during the preparation 
of this manuscript it was shown that the symmetry im¬ 
proved 2PIEA solves problems with infrared divergences 
of the standard model effective potential due to massless 
Goldstone bosons [3l| , though that study was carried out 
without the gauge sector. It also shows that the symme¬ 
try improved 2PIEA performs better than an ad hoc re- 


I I , ^ -•- -a 

a. — ^ig,_ <-a a-4 ^-a 4 

^ q = 0 
* 

^ --a H -^ ^-a H - A V-a 

I ) ^ ^ ^ I r ^ ^ 

1q = 0 q = 0 ^ 

I * >1 

1 q = 0 

FIG. VII.3. One loop contribution to the Higgs self-energy. 
The tadpole and sunset graphs are from the Higgs self-energy 
Eg, while the four remaining terms come from vertex cor¬ 
rections via the Ward identity (IIII.20II . The momentum in¬ 
coming from the lower Goldstone leg is zero, and the crossed 
vertex represents a factor of v. 


summation scheme proposed in the prior literature. This 
is heartening, though not wholly surprising due to the 
inherent self-consistency of nPIEA, a topic we plan to 
discuss in a forthcoming publication. 

However, the development of a first principles non- 
perturbative kinetic theory for the gauge theories of real 
physical interest requires the use of n-particle irreducible 
effective actions with n > 3. We have taken a step in this 
direction by extending the symmetry improvement for¬ 
malism to the 3PIEA for a scalar field theory with spon¬ 
taneous breaking of a global O (N) symmetry. We found 
that an extra Ward identity involving the vertex function 
must be imposed. Since the constraints are singular this 
required a careful consideration of the variational proce¬ 
dure, namely one must be careful to impose constraints 
in a way that satisfies d’Alembert’s principle. Once this 
is done the theory can be renormalized in a more or less 
standard way, though the counter-terms differ in value 
from the unimproved case. We derived finite equations 
of motion and counter-terms for the Hartree-Fock trun¬ 
cation, two loop truncation, and three loop truncation of 
the effective action. 

We found several important qualitative results. First, 
symmetry improvement breaks the equivalence hierarchy 
of nPIEA. Second, the numerical solution of the Hartree- 
Fock truncation gave mixed results: Goldstone’s theorem 
was satisfied, but the order of the phase transition was in¬ 
correctly predicted to be weakly first order (though there 
was still a large quantitative improvement over the unim¬ 
proved 2PIEA case). Third, the two loop truncation in¬ 
correctly predicts the Higgs decay width as a consequence 
of the optical theorem, though the three loop truncation 
gives the correct value, at least to 0{h). These results 
could be considered strong circumstantial evidence that 
one should not apply symmetry improvement to nPIEA 
at a truncation to less than n loops. One could test this 
conjecture further by, for example, computing the sym¬ 
metry improved 4PIEA. We predict that unsatisfactory 
results of some kind will be found for any truncation of 
this to < 4 loops. 

Our renormalization of the theory at two and three 
loops was performed in vacuum. The only finite temper¬ 
ature computation performed here was for the Hartree- 
Fock approximation. The extension of the two or three 
loop truncations to finite temperature, or an extension 
to non-equilibrium situations, will require a much heav¬ 
ier numerical effort than what we have attempted. It 
would also be interesting to compare the self-consistent 
Higgs decay rate in the symmetry improved 2PI and 3PI 
formalisms. We leave these investigations to future work. 
Similarly, we presented analytical results for the renor¬ 
malization of the three loop truncation only in 1 -|- 2 di¬ 
mensions, since the renormalization was not analytically 
tractable in 1 -|- 3 dimensions. This is also left to fu¬ 
ture work. The general renormalization theory presented 
here, based on counter-terms, is difficult to use in prac¬ 
tice. It will be interesting to see if symmetry improve¬ 
ment could work along with the counter-term-free func- 
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tional renormalization group approach (3^ . Such an ap¬ 
proach may not be easier to set up in the first place, but 
once developed would likely be easier to extend to higher 
loop order and n than the current method. Of course it 
will be interesting to see if this work can be extended to 
gauge symmetries and, eventually, the standard model of 
particle physics. If successful, such an effort could serve 
to open a new window to the non-perturbative physics 
of these theories in high temperature, high density and 
strong coupling regimes. 
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Appendix A: The d’Alembert formalism 

The assumption that Sf/ 6 W 2 can be consistently 
taken to be constant requires explanation. Constrained 
Lagrangian problems are generally under-specified un¬ 
less one invokes some principle like d’Alembert’s prin¬ 
ciple (that the constraint forces are “ideal”, i.e. they do 
no work on the system) to specify the constraint forces. 
Note that it while it is usually stated that enforcing con¬ 
straints through Lagrange multipliers is equivalent to ap¬ 
plying d’Alembert’s principle, this is no longer automat¬ 
ically the case if the constraints involve a singular limit 
as happens in the field theory case. This leads to a real 
ambiguity in the procedure which requires the analyst to 
input physical information to resolve it. In the case of 
mechanical systems the analyst is expected to be able to 
furnish the correct form of the constraints by inspection 
of the system. However, the interpretation of “work” and 
“constraint force” in the field theory case is subtle and 
the appropriate generalization is not obvious. Here we 
argue, by way of a simple mechanical analogy, that the 
procedure which leads to the maximum simplification of 
the equations of motion is the correct field theory ana¬ 
logue of d’Alembert’s principle in mechanics. 

d’Alembert’s principle is empirically verifiable for a 
given mechanical system, but for us it forms part of the 
definition of our approximation scheme, which we refer 
to as the “d’Alembert formalism.” The result of Sec¬ 
tion m was a set of unambiguous /-independent equa¬ 
tions of motion and constraint at some fixed order of the 
loop expansion, say Tloops. The use of any other limit¬ 
ing procedure requires the analyst to specify a spacetime 
function’s worth of data ahead of time, representing the 
“work” that the constraint forces do. The resulting equa¬ 
tions of motion represent a different formulation of the 
system and will have a different solution depending on 
the choice of “work” function. 

Imagine that we are competing against another analyst 
to find the most accurate solution for a particular system. 


It is possible that a competing smart analyst could choose 
a work function that results in a more accurate solution 
than ours, also working at Tloops. However, we could 
beat the other analyst by working in the d’Alembert for¬ 
malism but at higher loop order. We conjecture that the 
optimum choice of work function (in the sense of guar¬ 
anteeing the optimum accuracy of the resulting solution 
of the Lloop equations) is merely a clever repackaging of 
information contained in > Tloop corrections. (We have 
no proof of this conjecture. Indeed it is hard to see if any 
alternative to the d’Alembert formalism is practicable.) 
Thus we choose the d’Alembert formalism, which has the 
virtue of being a definite procedure requiring little clev¬ 
erness from the analyst, at the cost of potentially having 
a sub-optimal accuracy for a given loop order. 

To illustrate the connection with a mechanics problem 
consider a classical particle in 2D constrained to x'^+y^ = 
. The motion is uniform circular: 

The action is 

S = J Ldt-Xf [w], (A. 2 ) 

L = + y^) , (A.3) 

where the constraint is w (t) = x (t)^ -I- y (t)^ — = 0 

and / [w] = 0 if w (t) = 0. The equations of motion are 


mx (t) =-2Xj^^x (t), (A.4) 

my {t) =- 2 X-^^y {t). (A.5) 


In this mechanics problem we could set / [ui] = 
/ w (t) dt and carry through the problem in the standard 
way without any complications. But to mimic the field 
theory case, where a limiting procedure is required, we 
take f[w]= J w (t) dt. In this case 

^ . = 2w(t) ^ 0 as w ^ 0. (A-6) 

ow [t) 

This requires A —>■ oo such that X 6 f/ 6 w approaches a 
finite limit. Importantly, it must approach a t inde¬ 
pendent limit, otherwise an unspecified function of time 
enters the equations of motion: mx (t) = —k (t) x (t ), 
etc. This limit can be achieved by restricting the class 
of variations considered. Let x{t) = r (t) cos 6 (t) and 
y{t) = r (t) sm9 {t), where 6 r (t) =r{t) — r parametrises 
deviations from the constraint surface. Then w (t) = 
2r6r (t) -I- O (Jr^). We want w (t) = 0 which is obviously 
satisfied by Sr (t) = Sr. 

We are arguing that we only consider variations of 
this restricted form. The variations along the constraint 
surface (i.e. variations of 9 (t)) are unrestricted as they 
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should be. Only variations orthogonal to the constraint 
surface are restricted. This is equivalent to d’Alembert’s 
principle. To see this we compute the second derivative 
of to obtain 9^ — k = -. When the constraint is en- 

r 

forced f = 0, hence k ^ 0 implies 0^0: the constraint 
forces are causing angular accelerations, doing work on 
the particle. At constant radius, the centripetal force 
only changes if the angular velocity changes. 

In the field theory case we have pil.561) . For any given 
value of V and Aq only one value of Ah satisfies the 
constraint, given by 


(a;, y) = [ V (x, y,z)v + Aq^ (x, y) , (A.7) 

J Z 

where the * denotes the constraint solution. This is a 
holonomic constraint: in principle we could substitute 
this into the effective action directly and not worry about 
Lagrange multipliers at all (this is very messy analyti¬ 
cally, though it may be numerically feasible). We sug¬ 
gest that one restrict variations of A'Jj^ to be of the form 
Ajj^* {x, y) + Sk, where Sk is a spacetime independent 


constant. This way we guarantee 

——-- = 2 W 2 (x, y) = -25k = const, (A.8) 

5VV2 {x,y) 

and all the desired simplifications go through. Variations 
of the other variables are unrestricted. Because the con¬ 
straint force B5f / 5 W 2 disappears from the Aq, V and 
Vm equations of motion the constraint “does no work” 
on these variables, and the other variables {v and Ah) 
are determined solely by the constraint equations. This 
seems a fitting field theory analogy for d’Alembert’s prin¬ 
ciple. 

Appendix B: Deriving counter-terms for two loop 
truncations 

In this section we derive the counter-terms required to 
renormalize the 3PIEA and equations of motion in the 
two loop truncation as discussed in Section IIV Al Sub¬ 
stituting the expressions for bare fields and parameters 
in terms of the renormalized fields and parameters ac¬ 
cording to (IIV.9|) - (IIV.13I) gives the renormalized effective 
action 


r(3) = 




m 


dmn 


A 


n 2 T 4 

V --- V 


ih 




-y(iV-l)Tr 


2 ' 4! 

ZZ/^df^d^ + m'‘ + Srrii + Z^ 


+ - (V - 1) Tr In (Z-^Z^^Ac^) + -Tr In {Z-^Z^^Ajj^) 


2 , A-bdAf 2 


v^]Ag 


-^Tr 

2 


ZZ/^dnd^ + m'^ + i5to? -I- Za 


3A -b (5Af -b 25Xf 
6 




-I p(3) 
+ ^ 3 I 


(B.l) 


which agrees with the non-graphical terms of Q equation 
(4.4) upon setting N = 2, dropping an irrelevant con¬ 
stant oc TrlnZ“^ and noting our different conventions 
(m^gre “ —and Ahere = 6 Apt). The 5\ terms can be 
derived by substituting Xb^BcV’b (A + <5A(^) Zv'^ 


and XB^fBafBb —>■ Z ^ (a -b i5Af) Zv^SaNSbN into the 
definition of 

The graph functional becomes 


with 


,(3) _ ^ _ fi 2 ^A + dA^ 3 


rr = $1 - 


3! 




Ah (x, y) [A/f (x, z) Ah (x, w) Vn (y, -2, w) 


-b (V - 1 ) Aq (x, z) Aq (x, w) V (y, z, w)] - ^2 + 0 (ft^) 


(B.2) 


= ^ [(^ + 1) ^ + (^ - 1) ^>^2 + ^SX^] Zi {N - 1) AgAg 

+ ^ [3A + SX^ + 2SX^] ZIAhAh + ^2 (A + 5X^) Zl {N - 1 ) AgA^ (B.3) 

4)2 = ^ (A — 1) ZyZ\VVA hAgAg + —ZyZ^VnVnAhAhAh- 


(B.4) 
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The 6 X 2 ^^ terms can be found from substitut¬ 
ing XsABaa^Ebb {X + SX^) Z^Z^A^a^bb and 

XsABabABba Z {X 8 X 2 ^ Z"^ Z'^AabAba intO $ 1 - 
The $1 terms correspond to the Hartree-Fock approxi¬ 
mation and agree with the remaining terms of equation 
(4.4) of Q. The remaining O terms in Tg^^ give 


the sunset diagrams on replacing V and Vat by the solu¬ 
tion of their equations of motion at O {pP ^, which give 
Vn = 3V = -Zp^ (A + 6 X)vx (x - y) {x - z). 
We find that Zy cancels on elimination of V and Vn- It 
also disappears from the Ward identity once V is elimi¬ 
nated and hence plays no role in the further development. 

Going to momentum space the final result is (up to an 
irrelevant constant) 


r(3) = 


-z 


_impSrriQ 2 A-|-JAq 4 


A 


4! 


) + y(^-l)Trln(Ac') + yTrln {Ap^) 


ih 

T 


_ (TV - 1) / ( -ZZ^k^ +wP + 5ml + Za \ ^ ) Ac (fc) 


A-k(5Af 


- y y ^-ZZaA:^ + rrP‘ + 6m\ + Za 


3A-k 5Af-k 25Af 2 


Ah {k) 


+ ^ [(A + 1) A + (A - 1) 5X^ + 25Af ] Zi (^ - 1) ^ Ac (fc) Ac (fc) 

+ ^ [3A + 5X^ + 2(5Af ] Zi / A^ {k) Ah (k) + ^2 (A + Zi (A - 1) / Ac (fc) A^^ (k) 

J k k 

S2r(A + <5A)u] zl{N-l)[ AH{k)AG{l)AG{k + l) 

Jkl 


+ 


+ [(A + (5A) u] Za f Ah {k) Ah (1) Ah {k + 1). 

Jkl 


(B.5) 


From this expression we derive the renormalized equa¬ 
tions of motion (IIV. 1411 - (IIV. 161) 

The divergent integrals Tg/h (IIV-171) and Ihg (p) 
(IIV.181) enter into the equations of motion. Ihg {p) can 
be rendered finite by a single subtraction 

iHG{p)=i'^+:kfPG{p), (B.6) 

where [lA^^ ( 9 )]^- Since we wrote the propagators 

with the physical masses explicit, it is crucial to also 
subtract a portion of the finite piece ifpG /h) so that 
the pole of the propagator is fixed at the physical mass of 
the Goldstone/Higgs propagator respectively. We make 
this subtraction separately so as to have a universal X^. 

The tadpole integrals Tg/h require two subtractions 
each since i [A^ (g)]^ '^%/h ( 9 ) logarithmically diver¬ 
gent. To that end we introduce 


E^(q) = -ih 


(A -I- 5X) V 


n 2 




[ iA>^ (£) lA^^ {qx£)- x^ 

Ji 


(B.7) 


which is asymptotically the same as X^g/h ( 9 )’ that, 
i [A^ (g)]^ ^G/H ('(') ~ ( 9 ) finite. For later con¬ 


venience we write 


[ t [A'^ (g)]" (<z) = h 

J a 


(A -I- SX) V 


Zic^ (B.8) 


Then 


Tg/h = T‘^-^ (m^B - 


(A -I- i5A) V 


Zic^+^/H. 


(B.9) 


where = /^iA^(g). Note that and are real 
and X^ is imaginary, so that all of the subtractions can 
be absorbed into real counter-terms. 

The counter-terms are found by eliminating 
and demanding that the divergences proportional to dif¬ 
ferent powers of and TqJh separately vanish. Fur¬ 
ther, we enforce A^^ (me) = 0 and App {mn) = 0 
and that the counter-terms are momentum indepen¬ 
dent. This gives eight equations for the seven constants 
Z, Za, , i5A)^, <5A^, (5Af and i5A, however one of them 
is redundant and a solution exists [l^ . 

We find nontrivial field strength renormalizations 

Z = Z^^ = |l -I- -y [X^g {mn) — Tjhg (’^g)] | , 

(B.IO) 
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FIG. C.l. A contribution to the six point kernel K-abcdef re¬ 
sulting from (left to right) a derangement, two stabilizers and 
another derangement. The dashed boxes surronnd the per¬ 
mutations (to aid visualization only). Only stabilizers lead to 
divergent loops. 


and a nonzero 

(5A = - A ± A |l -b ^ ("iff) - | 

(B.ll) 

(the two solutions arise because (5A only enters the equa¬ 
tions of motion in the quadratic combination (A + 5X) ). 
These counter-terms are normally trivial {Z = = 1 

and 5X = 0) for theory at two loops. However, due 
to the modification of the Higgs equation of motion, we 
require Z\ 1 in order to enforce IS.^ {mn) — 0 and 
this is then compensated by Z and 5X in order to re¬ 
cover the other renormalization conditions. The other 
counter-terms can be obtained for any regulator, but the 
expressions are bulky and unenlightening even for dimen¬ 
sional regularization in d = 4 —2e dimensions, so we leave 
their explicit forms in the supplemental Mathematica 
notebook. 


Appendix C: Auxiliary vertex and renormalization 
in 3 and 4 dimensions 

As described in Section IIV HI the renormalization of 
the three loop 3PIEA requires the definition of an auxil¬ 
iary vertex with the same asymptotic behaviour as 
the full self-consistent solution at large momentum. This 
auxiliary vertex can be found in terms of a six point ker¬ 
nel which obeys the integral equation (IIV.29|) illustrated 
in Figure [TVT^ . 

Solving l|IV.29l) by iteration generates an infinite num¬ 
ber of terms, one of which is illustrated in Figure 1C.II 
Fach contribution is in one-to-one correspondence with 
the sequence of permutations 7 ri 7 r 2 • • • • • • of the prop¬ 

agator lines (read from left to right in relation to the dia¬ 
gram). Now we divide the permutations into two classes: 
“stabilizers,” for which tt (a) = a, and “derangements,” 
for which tt (a) = 6 or c. 

Any sequence of permutations is of the form of an al¬ 
ternating sequence of runs of (possibly zero) stabilizers, 
separated by derangements. Consider a run of n sta¬ 
bilizers, • • • tTq ( 7 ri 7 r 2 • • • 7r„) tt}, • • •, where tt^ and TTb are 
derangements and tti through 7 r„ are all stabilizers. The 
case for n = 2 is shown in Figure ICTl Fach stabilizer cre¬ 
ates a logarithmically divergent loop on the bottom two 
lines ~ —AZ^. Derangements on the other hand, if they 



FIG. C.2. Solution for the auxiliary vertex in terms of the 

four point kernel which sums all iterated bubble insertions. 


create loops at all, create loops with > 2 propagators, and 
hence are convergent. Thus all divergences in ICabcdef can 
be removed by rendering a single primitive divergence fi¬ 
nite. Note that the whole series '' ’^a (Iir^i ’’’0 

where again TTa,b are derangements and {iTi} are stabiliz¬ 
ers, can be summed because the series is geometric. The 
result is that the six point kernel can be determined by 
an equation like (IIV. 291 ) , except that the sum over all per¬ 
mutations is replaced by a sum over derangements only, 
and the bare vertex W is replaced by a four point kernel 
^itcd A/ (I -I- XI^). Denoting this four point kernel by 
a square vertex we can finally write the solution for 


in Figure 1C.21 

This expression for can be dramatically sim¬ 

plified in 3 or I -|- 2 dimensions because 
is finite and the geometric sum in JCabcd con¬ 
verges. Indeed (Pi>P 2 ,P 3 ,Pi + P 2 - Ps) ~ 

A/ l-\- X/ (pi — >■ A as Pi,2,3,4 —>■ 00. Further, 

every loop integral in Figure IC .2 1 likewise converges, and 
every loop yields a factor of ^ Thus the domi¬ 

nant behaviour as p —> oo is just the tree level behaviour 
and we can eliminate the auxiliary vertex completely. 


However, in 4 or 1 -b 3 dimensions apparently can¬ 
not be simplified further. First K-^^bcd "lust be renormal¬ 
ized, then the bubble appearing in the nontrivial terms in 
Figure (or the equivalent integral equation) must be 
renormalized, then the resulting series must be summed 
(or the equivalent integral equation solved), noting that 
on the basis of power counting every term is apparently 
equally important. On this basis we expect that no com¬ 
pact analytic expression for or even its asymptotic 
behavior, exists and that the renormalization must be 
accomplished as part of the self-consistent numerical so¬ 
lution of the full equations of motion. 

This style of argument can be quickly generalized to 
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many other theories, such as gauge theories, where the di¬ 
agrammatic expansion has a similar combinatorial struc¬ 
ture to scalar O {N) theory, showing up the well known 
problem of the renormalization of nPIEA for n > 3 in 
four dimensions. The discussion here certainly does not 
solve this problem, which remains open, to our knowl¬ 
edge, though we hope this discussion may be helpful. 


Appendix D: Deriving counter-terms for three loop 
truncations 

In this section we work in 1-1-2 dimensions as discussed 
in Section llVBI The effective action is as in Appendix [B] 
(before eliminating V and Vn) except we introduce a new 
counter-term 6X —^ SXc for the second term in (IB. 21) and 
add the three loop diagrams 


$3 = zt^zi 


- 1) (y)^ ^ + {N - 1 ) — (v) Ah^h (Ag)"^ 




3! 


4! 


$4 = 


ih^ (A -b SX) ^2 ^5 


24 


7^ 7^ 


2 (N — 1) VVn AqAq + — l) VVAh (Ag)^ + ^VnVn (Ah)^ 


+ 2 " {N-l)VV{AGfAHAH 


ih? (A -b i5A)^ 4 
=-lii- 


{[{N - 1) AgAg + AnAnf + 2 (TV - 1) (Ag)" + 2 (A^)^} . 


The equations of motion following from T*^^) are then 


(D.l) 

(D.2) 

(D.3) 


Ag,^ = — ( ZZ/^dfj,d^ -b -b Smf + Za 


2 , x^2 , V. A-b(5A^^_, 


6 


- ^ [(A + 1) A + (A - 1) SX^ + 2SX^] ZITg “ ^ + 5X^) ZITh 


— ihZyZ\ 


o (A -b 5Xc) Zy^v ^ ^ 


AhAgV 


+ h^Zt^Zl Vh {V) {Anf (Ag)' + (V) Ah Ah {AgY 


(A -b 5X) ^2 ^5 

--- Zj\rZj 


■‘V^A 


VVn (Ah)^ Ag -b (A -b 1 ) VV Ah {AgY + (AgY Ah Ah 


h^x + sxY^, 


18 

for the Goldstone propagator, 




(A-bl)(AG) -bAffA^^AG 


(D.4) 


W _ (^ + V 

V - ^ 


+ ihZYzl [Vn (V) (AhY Ag + {V) Ah (Ag) 
ih (A -b 5X) 2 


-b 


6 


Vn (AhY + (A -b 1 ) a (Ag)^ -b WAgAh 


(D.5) 


for the Higgs-Goldstone-Goldstone vertex. 


Vn — — (A -b SXg) vZy 


-1 


+ ihZy Z^ 


(A-1) (y)"(AG)' + (VAr)"(A^4)' 


-b 


ih{X + 5X) 


{N — l)V Ag Ag -b 3VAr {AhY 


(D. 6 ) 


for the triple Higgs vertex, and finally 

0 = Ag^ {p = 0)'!;, 


(D.7) 
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0 = ZyZ^V (p, -p, 0) + Ac' (p) - (p), (D.8) 

for the Ward identities. 

Note that the only divergent integrals in these equa¬ 
tions are the linearly divergent tadpole integrals Tq/h 
and the logarithmically divergent BBALL integrals (last 
line of (ID. 41) 1. By power counting with reference to Fig- 
ure llV.3l one finds that the third, fourth, and fifth lines of 


(ID. 41) produce finite self-energy contributions with lead¬ 
ing asymptotics ^ p“^, p~'^, and p“^ respectively. We 
can separate finite and divergent parts of Ag' as 

- [E°c (p) - {me)] - (p), (D.9) 

where 


-E^ (p) = -^ (iv +1) A (Tg - rn - (Th - rn 


and 


— ih 


o (A -I- SXc) Zy V ^ ^ 


AhAgV + 


Vn (V)'^ (Ah)^ (Ag)'^ + (V)'^ Ah Ah (Ag) 




-h 


r (A -H SX) Z\ 


2 r 


VVn (Ah)^ Ag -I- (A -I-1) VV Ah (Ag)^ + 3DD (Ag)^ Ah Ah 


2 ry^ 


(A -h (5A)^ Zi 
18 


{N + 1) (Ag)" -f AhAhAg - (A + 2) 


(D.IO) 


-E^ (p) = -E^ (tug) - (^{ZZa - 1) d^d>^ + 5m\ + + (Za - 1) 


_^(A+l)Ar'^-J [(A-l)<5A^ + 2(5Af]rG-^ [(A-f l)A + (A-l)a^ + 2af] (^i - l) Tg 
- JaT'^ - ^SX^Th - J (a + <5A^) {Zl - 1) Th 


- ih (Z^Zi - 1) 


^ (A -I- 6Xg) Zy^v ^ ^ 


AhAgV 

r-.\4 


+ ft" {Zt^Zl - 1) Vh {vy {Anf (Ag)' + (D) A^A^ (Ag) 


18 


VVn (A_f/)" Ag -I- (A -|- 1) VV Ah (Ag)" + 8VV (Ag)' Ah Ah 


(D.ll) 


are the finite and divergent parts respec¬ 
tively and we introduced the bball integral 
S'" = X^pA'"(g)A'"(p)A'"(p-h 9 ). In this split 

we have already assumed that (A-I-^Ag)^^' and 
(A -I- <5A) Z^ are hnite, which will turn out to be the 
case. Renormalization requires Eg (p) = 0. Note the 
explicit subtraction of E^ (wg) in order to fulfill the 
mass shell condition. Doing the same now for A^' we 
find the pole condition 


0 = ZyZ^V {niH, -ruH, 0) u -I- m]j - uig - Eg {uih) , 

(D.12) 


which requires 


ZvZn 


m% + ^G (n^g) ~ ntg 

V (mH,-mH,0)v 


(D.13) 


which is finite. We take for our other renormalization 
conditions the separate vanishing of kinematically inde¬ 
pendent divergences, implying 


ZZa = 1, (D.14) 

Sm! = -S°G (^g) - ^ (A + 2) AX'" + (A + 2) 

(D.15) 


(D.16) 
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II 

(D.17) 

Zy = , 

(D.23) 

0 = {N-1) SX^ + 2(5Af 


a 

1 

II 

(D.24) 

+ [{N + 1)X + {N- 1) dX^ + 2SXi] 

(zi-i), 

Z = 

(D.25) 


(D.18) 

-< 

1 

II 

(D.26) 

0 = 6 X 2 J- (a J- 6 X 2 ) (.^A ~ 1) ■ 

We also choose the conditions 

(D.19) 

SXc = - 1) A. 

(D.27) 


(A + 5A) Zl = A, (D.20) 

(A + ^c) = A, (D.21) 


to recover the tree level asymptotics for V and Vn- These 
conditions give a closed system of nine equations for the 
nine quantities Z, Z^ , Zy, (5Af, 6 X 2 ^^, SX, and 

SXc- 

These conditions determine 
a^ = (5Af = -^-^%^ = («2-1 )a, (D.22) 


Note that if k = 1 all of the counter-terms ex¬ 
cept 6 mi vanish. This is a manifestation of the super- 
renormalizability of (j)^ theory in 1 -|- 2 dimensions. The 
non-zero, indeed finite, values of all of the other counter¬ 
terms are not required to UV-renormalize the theory, but 
only to maintain the pole condition for the Higgs propa¬ 
gator despite the vertex Ward identity. 
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